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(54) Active acoustic noise reduction system 

(57) An active acoustic noise reduction system 
which comprises a single input transducer and an output 
actuator that are physically located next to each other 
in the same location, In one embodiment, the input 
transducer and the output actuator are a hybrid repre- 
sented by a single element. The active noise reduction 
system is located as close as possible to the noise 
source and functions to generate an antinoise cancella- 
tion sound wave with minimum delay and opposite 
phase with respect to the noise source. The noise re- 



duction system also comprises a non linearity correction 
circuit, a delayed cancellation circuit and variable gain 
amplifier. The system provides user control of the quiet 
zones generated by the system by varying the gain of 
the variable gain amplifier. The system provides a user 
with the ability In one embodiment, an echo canceler is 
utilized to remove echoes fed back from the output ac- 
tuator. In another embodiment, an input decoder is used 
instead of an echo canceler to remove feedback picked 
up from the output actuator 
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Description 

FIELD OF THE INVENTION 

5 [0001] The present invention relates generally to noise reduction systems and more particularly relates to acoustic 
noise reduction systems adapted to reduce the noise at a point relatively far from the noise source. 

BACKGROUND OF THE INVENTION 

10 [0002] Digital adaptive reduction of noise in the time domain is typically performed by sampling the analog output of 
a microphone that is appropriately positioned to sense the input noise. The sampled analog noise is then converted 
to digital format via an A/D converter, passed through an adaptive digital filter and then converted back to analog via 
a D/A converter before being output to a speaker. The analog output of a microphone is utilized as the input to the 
internal adaptive algorithm within the prior art noise reduction system. 

15 [0003] Since the group delay of such a system as described above is relatively large, a typical method of noise 
reduction employed in the prior art using such a system is to reduce the effect of periodic noise sources as opposed 
to arbitrary noise sources. 

[0004] Typical prior art noise reduction systems utilize an adaptive digital filter in the main data path to reduce the 
effect of the noise source. This causes the group delay of the system to be relatively large. 
20 [0005] For example, U.S. Patent No. 5,553,154, issued to Tamamura etal., discloses an adaptive filter that receives 
pulses that are synchronized with the period of the noise. The interval of the input noise determines the length of the 
taps in the delay line of the adaptive filter. 

[0006] U.S. Patent No. 5,613,009, issued to Miyazaki et al., discloses a vibration control system that a reference 
signal from the vibration source and an error signal from an object are input to the adaptive filter. Feedback means 
25 generates a feedback control signal in combination with a feedfoward signal, forms a drive signal for a vibrating means. 
[0007] U.S. Patent No. 5,627,896, issued to Southward et al., discloses an active system for controlling noise. The 
system operates by limiting the output gain from a gradient descent algorithm. 

[0008] Frequency domain analysis of the input data in the main data path of the system is another common technique 
utilized by prior art noise reduction systems. For example, U.S. Patent No. 5,271 ,062, issued to Sugita et al., discloses 
30 a noise reduction system that generates a signal having the same frequency but inverted phase relative to the input 
noise. 

[0009] U.S. Patent No. 5,347,586, issued to Hill et ai., discloses a noise control system comprising a reference 
microphone for generating a reference signal that is correlated with the noise emanating a noise source, a plurality of 
loudspeakers, a plurality of error microphones for generating a plurality of error signals and an adaptive controller. 
35 [0010] U.S. Patent No. 5,365,594, issued to Ross et al., discloses a vibration control system that utilizes a vibration 
input signal derived from a sensor sampling the vibrations to generate an output representative of the interference free 
vibration of a primary source. 

[0011] U.S. Patent No. 5,519,637, issued to Mthur et al., discloses an active structural acoustic control method for 
reducing the sound emitted through a structure. The method uses an array of transducers placed in a far field structure 
40 and an array of actuators mounted or embedded in that structure. Each is controlled by the system controller and uses 
reference signals derived from the noise source. 

[0012] Alternatively, another method of noise cancellation used in prior art systems places the microphone as close 
to the noise source as possible and the loudspeaker relatively far from the microphone so as to create a delay equal 
to the time for the noise to travel from the microphone to the speaker. This delay is intentionally created in order to 
45 match the internal signal processing time of the noise reduction system. The propagation time for the noise is configured 
to roughly match and compensate for the signal propagation time within the noise reduction system. This noise reduc- 
tion method is particularly useful for cancellation of noise in a duct such as an air conditioning duct. The internal signal 
processing is performed during the time that it takes for the sound waves to travel from the microphone to the loud- 
speaker. 

so [0013] Another prior art noise reduction technique related to the one just described, is to place the speaker close 
the noise source rather than far away from it, place a second microphone in the desired quiet zone and to adapt a 
digital filter utilizing the second microphone output. However, this method is useful for canceling repetitive noise only. 
[0014] The implementation and calibration of noise reduction systems made according to the prior art techniques 
described above are typically very difficult and correspondingly costly due to the need to supply a plurality of different 

55 elements, e.g., two microphones and a loudspeaker, and place these plurality of elements in their appropriate places. 
[001 5] Another disadvantage of prior art noise reduction systems is that if the canceled noise is not periodic in nature, 
then the physical dimensions of the cancellation system become relatively large. When the loudspeaker is located far 
away from the noise source , then the noise cancellation is localized and limited to a relatively small area. It would thus 
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be desirable to have a noise reduction system in which the physical design and calibration limitations of prior art noise 
reduction systems are avoided. 

[0016] U.S. Patent No. 5,410,607, issued to Mason et al., discloses a noise cancellation system for reducing noise 
radiated from a complex vibrating surface. The system includes a motion sensor, a controller having fixed transfer 
5 function and operative to generate an antinoise signal, and an acoustic driver operative to generate an acoustic anti- 
noise field that is substantially 1 80 degrees out of phase with the original noise field. The antinoise field reduces the 
original noise field using destructive interference. 

[0017] U.S. Patent No. 5,618,010, issued to Pla et al., discloses an active noise reduction system that utilizes a noise 
radiating panel to generate a noise that cancels the noise source. A piezoceramic actuator is connected to the panel. 
10 A sensing device such as an accelerometer is used in generating an opposite noise signal with which to cancel the 
noise source in the vicinity of the panel. 

[0018] U.S. Patent No. 5,343,713, issued to Okabe et al., discloses active noise control for three dimensional space 
including noise detecting means, sound generating means, noise control effect detecting means and a control device 
consisting of a phase of a wave signal fed from the noise detecting means wherein a plurality of noise control effect 
15 detecting means and sound generating means are disposed at a distance shorter than a wavelength of the noise to 
be controlled. 

SUMMARY OF THE INVENTION 

20 [0019] The present invention is an active acoustic noise reduction system which comprises an input transducer and 
an output actuator that are physically located next to each other in the same location. In one embodiment, the input 
transducer and the output actuator are a hybrid represented by a single element. The active noise reduction system 
is located as close as possible to the noise source and functions to generate the cancellation sound wave with minimum 
delay with respect to the noise source. 

25 [0020] The active noise reduction system, located very close to the noise source, functions to generate synthetic 
sound waves having a phase opposite that the noise. Both the noise source and the active noise reduction system 
might be situated within an enclosure or may be situated external to an enclosure. In one embodiment of the invention, 
the noise sound wave and the cancellation sound wave spread almost from the same point producing a high amount 
of noise cancellation. The output power of the cancellation signal is chosen so as to achieve maximum cancellation of 

30 the noise sound. 

[0021] Another embodiment of the invention is applicable when it is possible to place the noise reduction system 
very close to the noise source but the noise source body is much larger then the noise cancellation system. In the case 
when the noise body generates noise having the same phase in a direction towards the noise reduction system, good 
cancellation is achieved in the far field. Another benefit is that the noise cancellation is based on detecting the noise 

35 at only one point on the noise source body. 

[0022] Another embodiment of the invention is applicable when it is not possible to get close to the noise sound 
source or when the noise does not emanate from a single point or when the noise source body is relatively large 
compared to the system. In this case, a plurality of noise reduction systems are placed side by side, i.e., in an array 
configuration, to produce a 'wall of silence' with each noise reduction system generating cancellation sound waves in 

40 accordance with the noise detected at that particular point. Each element of the noise reduction array in this case 
operate independently of each other as opposed to each element being connected to some central processing unit. 
[0023] In another embodiment of the invention, the noise cancellation system is utilized as a speaker or sound gen- 
erating device for multimedia sources. In this mode, the noise cancellation system not only serves to perform noise 
cancellation, but also serves to generate meaningful sound. In other words, the system can function to both remove 

45 unwanted noise and intentionally insert sound such as music, thus replacing noise with music. This functions to reduce 
the annoyance of background noise by adding pleasant sounding music, for example. Another application to that em- 
bodiment is replacing the disturbing noise with significant sound in order to imply a known brand sound to equipment. 
[0024] The acoustic cancellation method of the present invention is based on the behavior of acoustic beam patterns 
in air. Cancellation of the noise is achieved in an area far from the noise source while in an area relatively close to the 

50 noise source there may be pockets of noise that exist. The length of the quiet zone, as measured from the noise source, 
is determined by the power of the cancellation signal generated and output by the system. Since the output acoustic 
beam pattern is dependent on the characteristics of the output actuator and on the main cancellation frequency that 
is used, the type of output actuator or the angle between a plurality of actuators may need to be varied in order to 
achieve optimum results for different noise frequencies. The noise reduction method of the present invention is capable 

55 of achieving effective cancellation of the noise when the surface of the noise source is complex given that the distance 
from the noise source to the point of cancellation is bigger then the length of the noise source itself. 
[0025] In addition to sensing sound from the noise source, the system also detects the sound from the output actuator. 
The portion of the input signal that is due to the output actuator is removed using an echo cancellation technique. If 
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the output and input transducers are acoustically separate elements and there exists acoustic delayed feedback in the 
system, then using an echo cancellation system is preferred. Another advantage of the echo cancellation system is 
the elimination of feedback sound emanating from walls, furniture, etc. and sensed by the input transducer. If there is 
no delayed time feedback from the output transducer to the input transducer and a directional input transducer is used, 
5 then a computation may be performed on the input signal, instead of using an echo cancellation system, to discern 
the actual noise signal from the input signal. 

[0026] In addition, the cancellation signal generated by the output actuator may be reflected from the noise source 
itself thus adding to the amount of noise present. In order to eliminate this type of noise, a delayed cancellation signal 
is generated by the system. 

10 [0027] The delay and phase shift applied to the cancellation signal is matched to the delay and phase shift associated 
with the reflection and feedback of the sound from the output actuator. 

[0028] There is therefore provided in accordance with the present invention an acoustic noise reduction system for 
reducing the effects of a noise source, comprising input transducer means for sensing the acoustic noise field generated 
by the noise source and for generating an input signal therefrom, output actuator means for generating an acoustic 

'5 output field that is effective to reduce the level of the acoustic noise field, correction means for adjusting the input signal 
generated by the input transducer to compensate for the non linear characteristics of the input transducer and output 
actuator, echo cancellation means for removing from the input signal a portion of the output of the output actuator 
means fed back through the input transducer means, the output of the echo cancellation means representing a signal 
corresponding to substantially the noise source by itself, antinoise means for generating an antinoise signal opposite 

20 in phase to the input signal, the output actuator means generating the acoustic output field from the antinoise signal 
and wherein the input transducer means is located in relatively close proximity to the output actuator means. 
[0029] There is also provided in accordance with the present invention an acoustic noise reduction system for re- 
ducing the effects of a noise source, comprising input transducer means for sensing the acoustic noise field generated 
by the noise source and for generating an input signal therefrom, output actuator means for generating an acoustic 

25 output field that is effective to reduce the level of the acoustic noise field, correction means for adjusting the input signal 
generated by the input transducer to compensate for the non linear characteristics of the input transducer, input de- 
coding means for removing extraneous signals from the input signal so as to generate a signal corresponding to sub- 
stantially the noise source alone, antinoise means for generating an antinoise signal opposite in phase to the input 
signal, the output actuator means generating the acoustic output field from the antinoise signal and wherein the input 

30 transducer means is located in relatively close proximity to the output actuator means. 

[0030] The correction means comprising storage means for storing a plurality of coefficients, coefficient processing 
means for dynamically updating the values of the plurality of coefficients stored in the storage means and means for 
generating a corrected input signal from the contents of the storage means and the input signal. 
[0031 ] The correction means comprising storage means for storing a plurality of coefficients, sigma generating means 

35 for outputting a signal corresponding to substantially the noise source only, coefficient processing means for dynami- 
cally updating the values of the plurality of coefficients stored in the storage means and means for generating a corrected 
input signal from the contents of the storage means and the input signal. 

[0032] The echo cancellation means comprises a digital filter having a delay line with a number of taps whose total 
delay time is equivalent to at least a system time delay of the noise reduction system, adaptation means for dynamically 
40 adjusting the coefficient values associated with each of the taps of the digital filter and summing means for adding the 
output of the digital filter with the output of the correction means. 

[0033] The antinoise means comprises a variable gain amplifier operative to generate an amplified signal 180 de- 
grees opposite in phase from the input signal and gain control means for dynamically controlling the gain of the variable 
gain amplifier. The gain control means is adapted to receive a manual input control signal from a user which determines 
45 the gain of the variable gain amplifier, the user able to vary the location of a quiet zone generated by the system by 
varying the input control signal. The input control signal is generated by the user remotely from the system and trans- 
mitted to the system via wireless communication means. 

[0034] The system further comprises a low pass filter operative to reduce oscillations present in the system derived 
from feedback of the acoustic output field to the input transducer. 
so [0035] Also, the system further comprises delay cancellation means for reducing the effect of echo signals caused 
by the antinoise means sensed by the input transducer. The delay cancellation means comprises a plurality of delay 
cancellation circuits wherein each delay cancellation circuit is operative to reduce the effect of the echo caused by 
previous delay cancellation circuits. 

[0036] Further, there is provided in accordance with the present invention a method for reducing the effects of a 
55 noise source, comprising the steps of sensing the acoustic noise field generated by the noise source and generating 
an input signal therefrom, generating an acoustic output field that is effective to reduce the level of the acoustic noise 
field, adjusting the input signal generated by an input transducer to compensate for the non linear characteristics of 
the input transducer and an output actuator, removing from the input signal a portion of the output of the output actuator 
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fed back through the input transducer, generating a signal corresponding to substantially the noise source by itself and 
generating an antinoise signal opposite in phase to the input signal, generating the acoustic output field from the 
antinoise signal. 

[0037] Also, there is provided in accordance with the present invention a method for reducing the effects of a noise 
5 source, comprising the steps of sensing the acoustic noise field generated by the noise source and for generating an 
input signal therefrom, generating an acoustic output field that is effective to reduce the level of the acoustic noise 
field, adjusting the input signal generated by an input transducer to compensate for the non linear characteristics of 
the input transducer, removing extraneous signals from the input signal so as to generate a signal corresponding to 
substantially the noise source alone and generating an antinoise signal opposite in phase to the input signal, the output 
10 actuator means generating the acoustic output field from the antinoise signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0038] The invention is herein described, by way of example only, with reference to the accompanying drawings, 
15 wherein: 

Fig. 1 is a schematic diagram illustrating the noise reduction system of the present invention applied to an example 
area having a noise source; 

Fig. 2 is a schematic diagram of a single noise reduction system applied to reduce the effect of a noise source 
20 showing the acoustic beam patterns generated; 

Fig. 3 is a schematic diagram of a plurality of noise reduction systems applied as an array of elements to reduce 
the effect of a noise source showing the acoustic beam patterns generated; 

Fig. 4 is a block diagram of a first embodiment of the noise reduction system of the present invention utilizing echo 
cancellation; 

25 Fig. 5 is a flow diagram illustrating the calibration method of the delay cancellation circuit; 

Fig. 6 is a block diagram illustrating the echo cancellation portion in more detail; 

Fig. 7 is a block diagram illustrating the input transducer and output actuator implemented as a hybrid combination 
in a single element; 

Fig. 8 is a flow diagram illustrating the gain control method utilized in the first and second embodiments; 
30 Fig. 9 is a flow diagram illustrating the first calibration method associated with the first embodiment; 

Fig. 10 is a flow diagram illustrating the second calibration method associated with the first embodiment; 

Fig. 1 1 is a block diagram of a second embodiment of the noise reduction system of the present invention utilizing 

a computational method to reduce echoes and oscillations; 

Fig. 12 is a flow diagram illustrating the first calibration method associated with the second embodiment of the 
35 present invention; 

Fig. 1 3 is a flow diagram illustrating the echo removal method of the present invention utilized in the second em- 
bodiment; and 

Fig. 14 is a flow diagram illustrating the second calibration method associated with the second embodiment of the 
present invention. 

40 

DETAILED DESCRIPTION OF THE INVENTION 

[0039] A schematic diagram illustrating the noise reduction system of the present invention applied to an example 
area having a noise source is shown in Figure 1. The noise reduction system, generally referenced 10, is preferably 

45 placed very close to a noise source 24. The distance X is the distance between the noise source and the system 10. 
[0040] The smaller the distance X between the noise source and the system, the better the noise cancellation 
achieved. The system 10 comprises an input transducer 30 such as a microphone and one or more output actuators 
such as loudspeakers. In the example system shown in Figure 1 , three output actuator loudspeakers 32 are shown. 
The orientation of the output actuators is such that the sound waves generated by the output actuators cancel the 

50 noise source sound waves. 

[0041 ] The noise source 24 is shown generating acoustic sound waves 40 that are sensed by the input transducer 
30 in system 1 0. The width of the noise source is denoted by the value W and its length is denoted by the value Y. The 
system 10 is located from the noise source at a distance X. The width of the system is denoted by the value Z. 
[0042] The noise source and the system are shown in a typical application such as a living room environment in a 

55 residence. The living room area 12 comprises typical furniture found in a living room, for example, two chairs 14, 20, 
a coffee table 22 and sofa 1 8. A person 1 6 is shown seated on the sofa and positioned within the effective quiet zone 
of the system. Although the noise source 24 and the system 10 are shown as separate entities, an alternative is to 
place both the system 10 and the noise source within a single enclosure (not shown). If the noise source is placed 
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inside the enclosure, the enclosure is regarded as the noise source. The size of the noise source 24 influences the 
cancellation of the noise source acoustic waves. If the dimension Y of the noise source is large relative to the half 
wavelength of the shortest interference noise signal, then a second noise reduction system should be installed on the 
other side of the noise source in order to achieve cancellation on that side. 

5 [0043] The mechanical structure of the output actuator used in the system has an effect on the quality of the noise 
reduction achieved. This is especially true in the case when the width W of the noise source is bigger then the width 
Z of the system which is related to the length of the output actuators. Preferably, the system 1 0 is a symmetric structure 
built from many small output actuator elements wherein each output actuator is oriented in a different direction. The 
output actuator elements may be spread over the length of the housing enclosing the system in order to provide cov- 

10 erage up to 1 80 degrees. Note that only one input transducer 30 (e.g., microphone) is connected to the noise reduction 
system 10. Each of the output actuators drive the same phase of noise cancellation wave. The number of output 
actuators can vary in accordance with the particular application and may be reduced to one. Alternatively, the input 
transducer and the output actuator can be combined in a hybrid input/output element that is positioned against the 
noise source. 

15 [0044] In a case of a hybrid transducer and the use of a plurality of output actuators, only the central actuator is used 
as a transducer. In the case of both separate input/output transducers and a hybrid transducer, the input transducer 
30 must be oriented in a direction towards the noise source 24. If a plurality of output actuators is utilized, then the 
centrally located transducer is oriented in a direction opposite the noise source. The output actuators other than the 
centrally oriented one must be positioned so as to achieve good noise cancellation, especially when a small number 

20 of output actuators is used. 

[0045] The proper position for the output actuators can be calibrated either manually or automatically at the time the 
system 1 0 is installed. Automatic calibration can be performed utilizing two motors (not shown), each oriented to handle 
a different axis of motion, attached to the output actuators. 

[0046] The total response time of the system, from the time the noise sound waves reach the input transducer to the 
25 time the noise cancellation signal is output by the output actuator, is preferably as short as possible. The effective 
length X is increased by an amount equivalent to the system delay time for the highest frequency component of the 
noise source wave. 

[0047] In general, the system 1 0 detects the noise source using the input sensor 30, amplifies the noise with inverse 
polarity and outputs it through the output actuator 32. Since the output actuator and the input transducer are located 
30 vary close to each other, the contribution of the noise signal to the amplitude of the input signal is much lower than 
that of the output transducer. 

[0048] In some configurations the noise cancellation signal generated by the output actuator may be reflected from 
the noise source as shown by the dashed line 26, thus increasing the overall level of noise in the system. To eliminate 
this type of reflected noise, a delayed cancellation signal is output by the system. 
35 [0049] The acoustic noise cancellation method of the present invention is based on the behavior of the acoustic 
beam pattern in the air of the noise and the signal output by the system. 

[0050] A schematic diagram of a single noise reduction system applied to reduce the effect of a noise source showing 
the acoustic beam patterns generated is shown in Figure 2. Good noise cancellation is achieved in an area far both 
from the noise source and the input transducer/output actuator while in areas close to the noise source there may be 

40 zones having higher levels of noise. The distance Q from the noise source represents an area of relative quiet. The 
length of Q is determined by the output power of the noise cancellation signal output by the system . Since the acoustic 
beam pattern is dependent on the characteristics of the output actuators and on the main frequency of the cancellation 
signal, it may be necessary to use different output actuators or to change the angle between them in order to achieve 
optimum noise cancellation for noise sources having different frequencies. For noise sources with complex surfaces, 

45 jt is possible to achieve good noise cancellation when the distance from the noise source to the point of cancellation 
is larger then the noise source itself. 

[0051] With reference to Figure 2, a plurality of output actuators 32 is shown generating an anti-noise signal in 
response to a noise source 50. The sound waves 52 of the anti-noise combine with the original noise waves at a point 
farfrom the noise source thus creating a quiet zone. When a plurality of output actuators is used, each of the actuators 
50 creates its own anti-noise field oriented at a specific angle. The acoustic field of all of the output actuators acting 
together combine with the noise source acoustic field to create a high intensity area, a high cancellation area or quiet 
zone and a low cancellation area. 

[0052] The low cancellation area 52 is due to the effect of noise emitted from points such as 53 that are far from the 
system 10. Better noise cancellation is achieved if the noise source body generates homogenous noise, i.e., all of the 
55 points 53 generate noise having the same phase. The distance of the quiet zone from the noise source is dependent 
on the energy content of the anti-noise. 

[0053] A schematic diagram of a plurality of noise reduction systems applied as an array of elements to reduce the 
effect of a noise source showing the acoustic beam patterns generated is shown in Figure 3. This scheme is utilized 
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when the surface of the noise source is large and complex or when it is difficult to place the noise reduction system 
close to the noise source. Noise cancellation is optimum when the noise generated by the noise source at a particular 
point on the body of the noise source has the same phase as the noise emitted from other points on the body of the 
noise that are oriented in the same direction. 

5 [0054] Three noise reduction systems 60, 66, 70 are utilized in this example to reduce the effects of the noise source 
58. The output actuator 62 of the system 60 functions to generate an anti-noise sound field 64. The output actuator 
68 of the system 66 functions to generate an anti-noise sound field 70. Similarly, the output actuator 72 of system 70 
functions to generate an anti-noise sound field 76. Thus, an array of independent output actuators 62, 68, 72 is used 
to create a quiet area at a distance from the array. The effectiveness of the virtual wall of quiet generated by the system 

10 is determined by numerous parameters. Such parameters include the distance between the individual noise reduction 
systems, the mechanical and electrical characteristics of the output actuators in relation to the main noise frequency, 
and the combined energy generated by the output actuators. 

[0055] A high level block diagram of a first embodiment of the noise reduction system of the present invention utilizing 
echo cancellation is shown in Figure 4. The first embodiment, generally referenced 80, of the system comprises means 

is for sensing the noise source, generating an anti-noise signal and outputting this anti-noise signal through one or more 
output actuators. In addition to sensing the noise source, the input transducer also senses the anti-noise signal output 
by the output actuator. As described previously, the amplitude of the anti-noise is larger than the amplitude of the noise 
signal itself. Echo cancellation is utilized to cancel the portion of the input signal associated with the signal output by 
the output actuator. If the input and output transducers are acoustically separate elements and there exists acoustic 

20 delayed feedback in the system, then using an echo cancellation system is preferable. If there is no acoustic delayed 
feedback from the output actuator to the input transducer, then a computing algorithm may be utilized to extract the 
noise signal from the total input signal. 

[0056] The input portion of the system comprises an input transducer 84, anti aliasing filter 88, amplifier 90 and N 
D converter 92. The input transducer may comprise a microphone, which is preferably directional and exhibits a very 

25 short delay. Use of a directional input transducer directed towards the noise source minimizes the sensitivity to acous- 
tical inputs other then the noise source. It is also desirable for the microphone to filter the acoustic input to maximize 
the sensitivity for the frequency range in use. In cases when the surface area of the noise source is large and close to 
half of the wave length of the disturbing noise, or when a signal to noise improvement is needed, then the input trans- 
ducer might be a summation of few identical input transducers within that area. 

30 [0057] Input transducers suitable for use with the present invention include, for example, electromagnetic based 
transducers, mechanical accelerators, electrical accelerators, piezoelectric and piezoceramic elements, vibration sen- 
sors, a capacitance microphone, a silicon microphone, and an optical microphone. 

[0058] The analog signal output by the input transducer is then passed through an anti aliasing filter 88 to the analog 
amplifier 90. The anti aliasing filter 88 is a low pass filter (LPF) that is constructed to exhibit minimum delay and designed 
35 according to the sampling frequency. The fixed gain analog amplifier 90 is adapted to reduce input transducer sensitivity 
to the minimum needed in order to eliminate the response of the system to sounds originating from sources other than 
the noise source. The signal output from the fixed gain analog amplifier 90 is then converted to digital via A/D converter 
92. The A/D converter is sampled at a rate that is as high as possible, with a resolution of at least 12 bits with 16 bits 
being preferred. 

40 [0059] The digital data output from the A/D converter is then adjusted to compensate for non linearities in the output 
transducer. Compensation for non linearities include multiplying the digital data by a coefficient stored in a look up 
table (LUT) 97 via the multiplier 93. The coefficients are calibrated dynamically during operation of the system. A 
coefficient processor 99 functions to calibrate the LUT coefficients based on the digital data output of the A/D converter 
and the output of the summer 94. 

45 [0060] The output of the multiplier 93 is input to an echo canceler 95 which functions to remove the echo reflected 
back from the output actuator and picked up by the input transducer. The cancellation signal generate by the echo 
canceler 95 is added to the output of the multiplier 93 via summer 94. The output of the summer is input to the equalizer 
101 that comprises a digital filter for correcting the frequency response gain and group delay of the analog elements 
in the system, including the output actuator and the input transducer. If the input and output transducers delay is large, 

so then the equalizer 1 01 may function also as a short time predictor to predict the output signal. The equalizer causes 
the input signal having different frequencies to be generated at the output transducer after a fixed time delay. The 
output of the equalizer 94 is input to a low pass filter 1 00 which limits the maximum frequency within the system. LPF 
100 comprises a digital finite impulse response (FIR) filter or infinite impulse response (IIR) filter having a low latency 
in the pass band. The use of the digital low pass filter 1 00 having a very short delay is optional but is useful to limit the 

55 band pass of the system to help avoid high frequency oscillations caused by the effect of feedback. 

[0061] The output of the digital low pass filter 100 is input to a variable gain digital amplifier 108 whose gain is 
controlled by a gain control circuit 1 06. The variable gain amplifier functions to perform an inversion of the noise signal 
to generate an anti-noise signal which functions to cancel the effects of the noise. Thus, the gain control circuit 106 
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sets the gain of the amplifier to a negative gain value. In addition to generating the basic anti-noise signal, the amplifier 
1 08 functions to prevent oscillations from occurring in the system. This is achieved by controlling the gain of the amplifier 
108. The gain of the amplifier can be adjusted relatively slowly and does not have to be performed in real time. 
[0062] A second input to the gain control circuit 106 is a manual gain control input that is provided by a user. A user 

5 can interact with the noise reduction system 1 0 by adjusting the gain of the amplifier 1 08. Using a suitable input device 
such as a remote control, a user positioned such as shown in Figure 1 , can control the location of the quiet zone to be 
any distance from the noise source by adjusting the gain of the amplifier. The gain adjustment capability of the noise 
reduction system is meant to achieve maximum cancellation at any arbitrary point within the cancellation zone. The 
gain may be preset or adjusted manually, such as by remote control. The gain control circuit 106 also comprises 

10 protection against the gain being increased too high so as to cause oscillations. 

[0063] The level of gain required to create a suitable quiet zone of high noise sound cancellation is dependent on 
several factors, e.g., the particular noise source, acoustics and dimensions of the area, the position of the user, etc. 
For known fixed sources of noise, e.g., air conditioners, vacuum cleaners, etc., the gain adjustment of amplifier 108 
can be performed once at the time of installation of the noise reduction system. The gain control method is described 

15 in more detail hereinbeiow. 

[0064] The output of the amplifier 108 is input to a summer 112 and a delayed cancellation circuit 110. The output 
of the delayed cancellation circuit 1 1 0 is added to the output of the amplifier 1 08 via summer 112. The function of the 
delayed cancellation circuit is to generate a cancellation signal in response to reflections of the portion of the output 
signal from the output actuator (Figure 1 at 26) produced by the noise cancellation system itself. The reflection is 

20 effected by the acoustical environment, the distance of the system from the noise source 24 and the frequency of the 
noise. The delayed cancellation circuit is based on a delay circuit in combination with a negative gain multiplier. In the 
alternative, the circuit can be based on a digital filter having the same impulse response as the reflected signal but 
with the opposite, i.e., negative, polarity. 

[0065] The noise reduction system 1 0 effectively considers the acoustic waves reflected from the noise source as a 
25 second mirrored noise source located at a distance equal to the distance between the noise source and the noise 
reduction system. The only difference being that the path of the acoustic sound waves from the second noise source 
is in the opposite direction from the sounds waves emitted from the noise source itself. 

[0066] The output of the delayed cancellation circuit is added to the output of the amplifier 1 08. Note that the signal 
generated to cancel the first reflection may itself be reflected off the noise source to effectively create a third mirror 
30 noise source and so on and so forth. This third noise source can subsequently be canceled using the same method 
of the delayed cancellation circuit 110 but with different parameters (not shown). 

[0067] Note that the delayed cancellation circuit is optional but using it, however, can aid in canceling acoustic sound 
waves fed back from the noise source body. Calibration of the gain and the delay of the impulse response of the delayed 
cancellation circuit 1 1 0 is performed once during installation of the noise reduction system 1 0. A high level flow diagram 

35 illustrating the calibration method of the delay cancellation circuit is shown in Figure 5. The calibration procedure 
comprises first turning off the noise source and the echo canceler within the system (step 290). Then, a single impulse 
is generated and output through the output actuator (step 292). The impulse is generated towards the noise source 
while the noise source itself is off. The impulse response is then measured, i.e., the input fed back from the noise 
source is measured (step 294). The noise reduction system 10 regards this source of noise as a second noise source 

40 having known characteristics. 

[0068] The measured impulse response is then sampled (step 296) and FIR filter coefficients are generated based 
on the sampled impulse response (step 298). The gain of the delayed cancellation circuit 110 is determined using the 
following equation. 



so [0069] Wherein GAIN DCC is the gain of the delay cancellation circuit, GAIN mp is the gain of the variable gain amplifier 
1 08, IN PUT is the maximum amplitude of the response measured at the output of the A/D converter 92 before the data 
is corrected for non linearities and OUTPUT is the output impulse transmitted value of the variable gain amplifier 1 08. 
Thus, the gain of the delayed cancellation circuit (DCC) is based on the negative or inverted gain of the amplifier 1 08. 
[0070] Note that the initial portion of the measured impulse response is due to the input transducer picking up the 

55 output of the output actuator without echoes, reflections or feedback. The actual response to the impulse function does 
not begin until a certain time period after the impulse is generated. This time period is proportional to the system time 
delay period and can be calculated since the start of the impulse function and the measured response are known. The 
portion of the impulse response prior to the system delay period is discarded and not used to generate the coefficients 
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for the digital FIR filter that can be utilized in the delayed cancellation circuit. The system delay time is defined as the 
time from impulse generation until the time that the amplitude of the first response at the output of the A/D converter 
92 is equal to 0.2 of the maximum amplitude response. For example, if it is assumed that the system delay time is 1 00 
u.s, then the first 100 us of the measured impulse response is discarded. The portion of the response 100 u.s and 
5 beyond is sampled and used to generate the filter coefficients. If the system includes cancellation circuitry for the third 
mirrored noise source, then the calibration is performed in steps. In the first step, the second mirrored source cancel- 
lation coefficient calibration is performed as described above. With this cancellation activated, calibration of the second 
step is performed by repeating the procedure described in Figure 5. The new coefficients are derived from the second 
feedback filter (not shown). 

10 [0071 ] The output of the digital summer 1 1 2 is input to the D/A converter 1 1 4 and to the echo canceler 95 portion of 
the system. The echo canceler 95 comprises a digital filter 96, an adaptation circuit 98 and a summer 94. The echo 
canceler 95 functions to remove reflections of sound and prevent oscillations caused by the output of the output actuator 
feeding back to the input transducer. The filter 96 is a digital FIR filter having a sufficient number of taps, i.e., delay, to 
cover the round trip delay through the system as well as feed back echoes from the acoustical environment. 

15 [0072] In order to reduce the computation burden of the FIR, interpolation and decimation methods, well known in 
the digital filter arts, are used (not shown) at the FIR output and input, respectively. The coefficients of the FIR filter 
correspond to the coefficients of the impulse response of the portion of the system from the output of the summer 112 
to the input of the summer 94 and are measured and setup at the time the system is started up. The measurements 
is performed by using the maximum length sequence technique (MLS), a technique well known in the art. The MLS 

20 technique is implemented by generating the MLS pseudo random sequence atthe output of the summer 112, measuring 
at the input to summer 94, and then low pass filtering the cross correlation of the two signals. The results of the MLS 
measurement are the filter 96 coefficients. 

[0073] Real time adaptation of the digital filter coefficients by the adaptation circuit 98 is performed using any suitable 
FIR adaptation algorithm well known in the art, such as the least mean squares technique. The adaptation is based in 

25 the error signal output by the summer 94. 

[0074] Since the input noise signal from the noise source 24, e.g., motor or airflow noise, is a highly correlated signal, 
the regular FIR adaptation is applicable only outside of the auto correlation area. For example, in a typical air condi- 
tioning system, the auto correlation signal is 300 microseconds. Assuming a system delay of 1 00 microseconds, there 
are 300 minus 100 or 200 microseconds that may not be calibrated using the least mean square technique. 

30 [0075] If the length of the FIR is adapted to be 5 milliseconds, then only the last 4.8 milliseconds may be adapted in 
real time using the least mean square (LMS) technique using a special variance explained bellow. Note that the LMS 
technique is well known in the art. The coefficients within the first 200 microseconds change relatively slowly with time 
and are effected mostly by temperature and humidity. Calibration of the coefficients within the first 200 microseconds 
is performed slowly in real time using the MLS technique described above with the exception that no generation of 

35 signal occurs. 

[0076] Instead, the noise data is used to perform cross correlation. Averaging a large number of test results yields 
the coefficients for the first 200 microseconds of the FIR filter 96. Note that when using the local closed loop control 
system to control the linearity of the output actuator as described in the linearity table, there is no need to adapt the 
first 200 microseconds in real time since the system response during that time is constant. 

40 [0077] The noise is cyclic, i.e., repetitive, and to enable acceptable adaptation, only one noise cycle is enabled within 
the length of the FIR filter. The full length of the FIR filer is divided into a plurality of short FIR filters. For example, if 
the noise cycle is 600 Hz, then the maximum allowed FIR length is 1/600 = 1.6 milliseconds. To cover the full 4.8 
milliseconds, three short FIR filters are used wherein each filter spans 1.6 milliseconds. Each FIR filter comprises it's 
own separate adaptation circuit. 

45 [0078] The following example is regarded to LMS adaptation algorithm, but might be implemented on any known 
adaptation algorithm. The block diagram shown in Figure 6 illustrates the difference between a conventional LMS 
adaptation and the LMS adaptation of the present invention. Three FIR filter units are used with three taps in each FIR 
filter. The data from the summer 112 is input to the shift register 320 that functions as a serial to parallel converter. The 
first three registers 1 , 2, and 3 of the shift register 320 are utilized by the first FIR 321 in accordance with the following. 

so The input noise that comes from multiplier 93 is multiplied by the constant u.1 via multiplier 326. The resulting product 
is multiplied with the contents of the first three registers of shift register 320 via multiplier 325. The resultant product 
is subsequently added to the current coefficient 323 via adder 324. 

[0079] It is important to note that the size of the coefficient 323 and the size of the FIR filter is 1/3 of the total length 
of the FIR filter. The remainder of the length is used as a shift register delay 322. The FIR portion 321 and the shift 
55 register delay portion 322 combine to form the complete FIR length. The summer 329 sums the output of the FIR with 
the input from multiplier 93 to generate first part of the result. The output of the summer is input to the summer 343. 
[0080] The second set of three registers 4, 5 and 6 are delayed in sync with delay 330 which generates a delay equal 
to the length of FIR 321 . The output of the delay 330 is input to FIR 331 . The coefficient 334 adaptation is performed 
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using multipliers 328 and 327 and adder 333 in a manner similar to that of the adaptation of coefficient 323 but utilizing 
constant u2 and registers 4, 5 and 6. The output of the FIR 331 is input to delay 332 to complete the length of the FIR. 
The output of delay 332 is summed with the input from multiplier 93 via adder 335 to yield the second part of the result 
which is input to summer 343. 

5 [0081] The third portion of the output of summer 343 is generated utilizing bits 7, 8 and 9 of shift register 320 after 
passing the data through delay 340 and FIR 341 . The coefficient 339 adaptation is calculated using u.3, the input signal 
from multiplier 93, the output of multipliers 337, 336 and the output of adder 338. Note that the length of delay 340 is 
equal to the length of FIR 321 plus the length of FIR 331 . Thus, FIR 341 functions to complete the total length of the 
FIR filter. 

10 [0082] The dynamically changing values u.1, u2 and u.3 function to determine how fast the adaptation algorithm 
performs. For example, if the adaptation algorithm runs too fast, the LMS will not yield correct results, i.e., it may not 
converge. A typical value for u.1 , u.2 and u.3 is in the range of 0.1 to 0.2 which was derived from experimentation. 
[0083] The output of the summer 94 is a digital signal wherein the feedback, i.e., echo, caused by the generated 
output signal has been removed. The digital filter 96 functions to generate a signal that is substantially the opposite of 

is the portion of the input signal that is fed back from the output signal. 

[0084] The response of the adaptive digital filter 96 is the same response of that of the system from the output 
actuator to the input transducer including the environmental acoustic effects. The digital filter 96 functions to generate 
the signal output by the output transducer but with opposite polarity, i.e., negative polarity and the correct phase shift 
and amplitude. The summer 94 adds the output of the digital filter 96 with the output of the multiplier 93, which is the 

20 digital input corrected for the non linearities of the input transducer, to eliminate the echo from the input signal. The 
addition operation performed by summer 94 may be performed digitally as described above or can be performed in 
analog. The summing can be performed in analog by utilizing a D/A converter to convert the output of the digital filter 
96 to analog and summing the result with the analog signal from the input transducer. In the case of the analog echo 
cancellation summation, it is not economical to use an analog multiplier for the linearity table and thus, the assumption 

25 is that the actuator is linear. In either case, the output of the echo canceler is input to the digital low pass filter 1 00 as 
described hereinabove. 

[0085] It is important to note that the time response of the digital filter 96 must be long enough to cancel all of the 
acoustical sound fed back to the input transducer from the output actuator. The adaptive echo canceler 95 adjusts the 
filter taps in response to changes in the acoustical environment, e.g., temperature, furniture placement, etc. 

30 [0086] The output of the summer 112 is input to the D/A converter 114 which functions to convert the digital signal 
into analog. The output of the D/A converter 114 is input to an analog low pass filter 116. The cutoff frequency of the 
low pass filter 116 is calculated according to the D/A sampling rate and is designed to be low latency. The output of 
the low pass filter 1 1 6 is input to the power amplifier 118. The power amplifier boosts the level of the output sufficiently 
to drive the output actuator 120. The output actuator can be any suitable device that transforms the output signal to 

35 sound waves. The output actuator can comprise a loudspeaker having a relatively short delay. The output actuator 
120 may comprise a single output actuator or a plurality of output actuators connected in parallel. 
[0087] In addition , the system of Figure 4 also comprises a controller (not shown) which functions to administer and 
control the configuration, operation, settings and all timing of the noise reduction system 10. 
[0088] It is noted that the sound emitted from the output actuator is detected by the input transducer at a much higher 

40 amplitude than the noise itself. The noise reduction system 1 0 functions to discern the much smaller noise signal from 
the much larger echo signal. 

[0089] However, the input transducer exhibits a non linear response to the amplitude of the sound waves incident 
on it. It turns out that the input/output response of the input transducer, such as a microphone comprising a membrane, 
is dependent on the total amount of power hitting the microphone. Thus, for the same amplitude noise signal, the 

45 linearity of the microphone is different for an input with low total power input as compared to an input with high total 
power. Thus, the high amplitude signal output from the output actuator effects a different relative output signal from 
the input transducer due to the noise. This is especially true when using a hybrid to perform both the input and the 
output transducer functions of the system. Non linearity at the output actuator is manifested when the output peak 
power to the transducer is large. When using separate input/output devices, the non linearity in the output device is 

so much more critical when the mechanical amplitude is large. The operation of the system is based on the correct dis- 
crimination of the noise source from the much larger echo signal. Thus, a correction of the echo signal amplitude as 
it appears in the input is required in order to achieve good cancellation of the echo signal. The linearity LUT performs 
this correction operation by multiplying each input value with the correction coefficient. 

[0090] The linearity table may be generated a priori and calibrated at the time of installation or adjusted dynamically 
55 on the fly. The linearity table and the dynamic adjustment of its contents are described in more detail hereinbelow. 
[0091 ] As an alternative to the linearity table, a local closed loop control system can be used to compensate for the 
linearity of the actuator (not shown). Such a technique is described, for example, in the proceedings of Active 97 in 
the article by Yoon-Sun Kim and Youngjin Park entitled "Non Linearity Compensation for Harmonic Distortion of Direct 
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Radiation Loudspeaker." 

[0092] A high level block diagram illustrating the Input transducer and output actuator implemented as a hybrid com- 
bination in a single element is shown in Figure 7. As an alternative to a separate input transducer and output actuator, 
a hybrid input/output (I/O) element or transducer 1 62 can be utilized that performs the functions of both elements. 

5 [0093] In this case, a hybrid circuit 1 60 is needed to interface the I/O element 1 62 to the power amplifier 1 64 and 
the anti aliasing filter 1 66. The hybrid 1 60 functions to transfer power from the power amplifier 1 64 with minimum losses 
to the I/O transducer element 162 and with maximum reduction to the anti aliasing filter 166. 
[0094] The gain control method implemented in the gain control circuit 106 (Figure 4) of the variable gain amplifier 
108 will now be described in more detail. A high level flow diagram illustrating the gain control method utilized in both 

10 the first and second embodiments is shown in Figure 8. This method operates in a loop to continuously search for 
oscillations in the system. In general, an FFT or DFT performed in the gain control circuit is used to map the mean 
amplitude of the frequency content in the detected input. When a new frequency is detected or when the total input 
power to the system increases, an immediate reduction of the gain of the amplifier 108 is performed. 
[0095] The first step is to check the total power of the output signal to see if it is over a predetermined maximum 

15 (step 130). If the power is over the permitted maximum then the gain is reduced until the total output power is below 
the maximum (step 132). Once the gain is within the permitted range, the current gain setting is stored as the initial 
gain value (step 134). A Fast Fourrier Transform (FFT) or DFT is then performed on the signal input to the amplifier 
(step 1 36). The FFT or D FT yields a map of the frequency content of the input signal. A plurality of samples of the input 
are taken and FFT or DFT analysis performed on each sample. Corresponding frequency elements are averaged over 

20 many samples (step 138). 

[0096] It is then checked for the significant presence of signal content at new frequencies (step 1 40). If signal content 
is found at new frequencies, these frequencies elements are tracked overtime (step 142) to determine whether they 
are oscillations (step 144). A new frequency that is persistent in time is considered a suspect oscillation. If there is a 
suspected oscillation, then the gain is reduced by a step amount (step 146). If the gain is reduced to a predetermined 

25 minimum (step 148) the frequencies that were suspect as oscillations are mapped as input noise, i.e., a frequency 
generated by the environment with a particular amplitude and frequency range and not as oscillations (step 1 52) . This 
is because, at such low values of gain, it would be highly improbable that the suspect frequencies were caused by 
oscillations. Signals within the frequency range are treated as environmental noise. After the environmental frequencies 
are mapped, the gain is restored to the gain that was initially saved during step 134. 

30 [0097] If there are no new frequency elements in step 140 or oscillations suspected in step 144, the method ends. 
If the gain in step 148 is not a minimum, than another group of samples is input (step 150). After sampling another 
group of input signals, it is then checked whether oscillations are still present after the gain has been reduced. The 
method keeps looping by checking for oscillations, reducing the gain if oscillations are found and generating additional 
input sample until either no oscillations are found or the gain reaches a minimum. 

35 [0098] As described previously, the linearity LUT 97 is used to correct the non linear behavior of the input/output 
transducers. Since typical input transducers are constructed from moving mechanical parts, the output response for 
the same input noise level is different for different levels of total input power incident on the input transducer. However 
the output actuator suffers from non linearities when the motion of the moving mechanical parts increase as a result 
of increasing output power. The difference in output response of the output actuator and input transducer is corrected 

40 utilizing the linearity LUT 97. 

[0099] The linearity LUT stored different coefficients for each amplitude input value range. The linearity LUT is divided 
into regions with each region having its own coefficient value. For example, assuming a 12 bit data word, the linearity 
LUT may be divided into 256 regions, thus only utilizing 8 bits of the input data. Upon startup of the noise reduction 
system 10, all of the coefficients in the linearity LUT have the same value. 

45 [0100] The present invention includes a first and a second method of calibrating the linearity LUT. Both methods are 
performed in real time during the operation of the noise reduction system. Both calibration methods utilize a linearity 
look up table (LUT) that holds coefficients used to adjust the input data output by the A/D converter. The adaptation 
or calibration of the coefficients of the linearity LUT serves to compensate for slow changes to the linearity of the system 
caused by temperature, humidity etc. The linearity LUT also functions to compensate for the mechanical non linearities 

50 of both the input transducer and the output actuator(s). Although not necessary to perform the invention, the input 
amplitude is divided into a plurality of regions wherein each region has a coefficient associated with it. 
[0101] The first calibration method utilizes the input noise signal itself to calibrate the coefficients. In the second 
method, a small amplitude calibration signal is injected into the system during operation and the results used to generate 
new coefficient. Each calibration method will now be described in more detail below. The high level flow diagrams 

55 describing the linearity table calibration methods refer to the calibration of one coefficient. The methods described are 
repeated in order to calibrate all the coefficients. 

[0102] A high level flow diagram illustrating the first calibration method associated with the first embodiment is shown 
in Figure 9. The first method of calibration utilizes the fact that the noise is physical and continuous. The controller in 
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the system tracks the relationship between values termed Table Input fTI), Table Output (TO) and Summer Output 
(SO) during operation of the system. TheTI values are measured at the output of the A/D converter 92 (Figure 4), the 
TO values are measured at the output of the multiplier 93 and the SO values are measured at the output of the summer 
94. The coefficient processor 99 functions to calculate new LUT coefficients based on the Tl, TO and SO values, 
[0103) The calibration of the LUT coefficients during operation of the system attempts to ignore the effects of the 
noise source. Note that the input noise source itself changes between two adjacent samples. Note also that the output 
of the summer SO represents the noise source since the echo canceler 95 effectively removes the echo signal from 
the input signal. 

[0104] The following calibration method is based on the system time delay which is defined as the time from TO to 
the D/A 114 plus the time from the D/A 114 to the A/D 92. Note that factored in this time is the acoustic medium, the 
analog transducers the circuit and the time from the A/D 92 to TO. The time attributed to the acoustic medium and the 
analog transducers and circuit is substantially equivalent to the time during which most of the energy of the impulse 
generated at D/A 114 is measured at the A/D converter 92. 

[01 05] The measurement of the impulse is performed using the impulse response measurement technique described 
previously in connection with the calibration of FIR filter 96. 

[0106] The calibration method initially measures two TO values, i.e., a value at time 'n-1' and another at time 'n\ 
Subsequently, after a system time delay, the controller then measures the two Tl and SO values that are the effect of 
the previously measured TO values, assuming the system time delay is known. In other words, each Tl and SO value 
is measured after the corresponding TO values has had a chance to propagate through the system, i.e., output by the 
output actuator, sensed at the input transducer, etc. Note that the SO value is the Tl value after compensation for non 
linearities and after the echo fed back from the output actuator to the input transducer is removed to yield a value that 
reflects the noise level only. 

[0107] With reference to Figure 9, the first step in the calibration method is to read a Table Output (TO) value at the 
output of the multiplier 93 (Figure 4), termed TO r>1 (step 21 0). After reading the TO n _., value, the system waits for the 
output of the multiplier to change before reading the next value which is termed TO n (step 214). During operation of 
the system, it is not unusual if the input data does not vary much from sample to sample. This is due to the fact that 
the sampling rate of the system is orders of magnitude higher than the frequencies making up the noise signal. In step 
214, the system waits more then one sample time, and the sample is taken when the difference between the level of 
TO^ and TO n reaches some minimum which is predetermined. This minimum determines the coefficient calibration 
accuracy. It is then determined whether TO n is within the same region of the LUT as TO n .., (step 216). This is checked 
in order to prevent two TO values being associated with different regions of the LUT. The calibration method calculates 
new coefficients for a single region of the LUT at a time. The calculations, thus, cannot span borders between regions. 
[01 08] Next, the system waits for the effect of the TO n _ 1 value (step 21 8) and TO n value (step 222) to appear at the 
output of the A/D converter and the outp ut of the summer . As described previously, the data output by the A/D converter 
is termed Table Input (Tl) data and the data output by the summer is termed Summer Output (SO). Once the effect of 
the TO n _., data appears at the output of the A/D converter and the summer, the Tl n _ 1 and SO n .., values are read (step 
220) . Similarly, once the effect of the TO n data appears at the output of the A/D converter and the summer, the Tl n and 
the SO n values are read N n (step 224). The steps of first waiting and then reading the Tl and SO values described 
herein can be implemented either sequentially or in parallel. 

[0109] The index to the LUT which determines which coefficient is presently under calculation is then generated 
using the Tl^ value read during step 220 (step 225). This index is used for the calibration process only and does not 
effect the main real time data path of the system. In the example noise reduction system shown in Figure 4, the LUT 
has less entries in it than the number of possible input values, e.g., 256 regions for 12 bits of input data, in order to 
reduce the size of the lookup table required. Alternatively, the LUT can be constructed to hold a coefficient value for 
each and every possible input data. 

[0110] Once the Tl and SO data have been read, the new coefficients are calculated using the following equation 
(step 226). 



[0111] The new coefficient C new is a function of the old coefficient C dd . The values TO, Tl and SO are used to 
generate an intermediate new coefficient from which C oid is subtracted. A portion of the delta is added to C o(d to perform 
the calibration. The constant K varies between 0 and 1 and is used to determine the speed with which the coefficients 
are permitted to change. Values of K closer to 0 cause the coefficients to change more slowly whereas values of K 
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closer to 1 cause the coefficients to change more quickly. 

[0112] The second calibration method associated with the first noise reduction embodiment of Figure 4 will now be 
described in more detail. A high level flow diagram illustrating the second calibration method associated with the first 
embodiment is shown in Figure 10. This second method of calibration is similar to that of the first method described in 

5 connection with Figure 9, with the difference being that rather than wait for the actual noise source to cause a change 
to the TO value, an artificial noise signal is injected into the data path to simulate a known change in the noise signal 
level. The second method of coefficient calibration models the noise as a pseudo DC level noise source during cali- 
bration. This is a reasonable assumption since the calibration period of approximately 10 ^ts is very short relative to 
the noise frequency. This pseudo DC level of the noise is used to point to a particular region in the linearity LUT 

10 [0113] With reference to Figures 4 and 1 0, the first step is to measure the Table Input (Tl) value denoted Tl n . 1 at the 
output of the A/D converter 92 (step 230). The index to the linearity LUT is then generated based on the Tl fv1 value 
just measured (step 232). The index determines which of the coefficients of the linearity LUT is to be calibrated during 
this particular invocation of the method. The next step is to read the Table Output (TO) value, denoted TO^, at the 
output of the multiplier 93 and the Summer Output (SO) value, denoted SO^, at the output of the summer 94 (step 

15 234). The TO value is generated by multiplying the output of the A/D converter with the output of the LUT The result 
of the multiplication is added to the output of the digital filter. Note that the TO and SO values are read immediately 
after the Tl value is read without waiting a system time delay as in the first method of Figure 9. 
[0114] A calibration signal is then injected at the output of the multiplier (step 236). The output of the multiplier which 
is denoted as the TO value is replaced with the calibration signal for a finite time period. The calibration signal, termed 

20 TO n , comprises the original output of the multiplier TO n ..j increased by a known delta amount. The system then waits 
one system delay time for the injected calibration signal to appear at the output of the A/D converter (step 238). After 
waiting one system time delay, the data at the output of the A/D converter is read and termed Tl n . In addition, the output 
of the summer 94, termed SO n , is also read (step 240). 

[01 15] After reading the Tl n value, the original output of the multiplier TO n _.| before the calibration signal was injected 
25 is restored (step 242). Then, based on the values Tl, TO and SO, the new coefficient can be calculated utilizing the 
following equation (step 244). 



30 



C = C old + K 



TO - TO.. 



iTI^SOj-iTI^-SO^) 



[0116] The new coefficient C new is a function of the old coefficient C oJd . The values TO, Tl and SO are used to 
35 generate an intermediate new coefficient from which C old is subtracted. A portion of the delta is added to C 0 , d to perform 
the calibration. The constant K varies between 0 and 1 and is used to determine the speed with which the coefficients 
are permitted to change. Values of K closer to 0 cause the coefficients to change more slowly whereas values of K 
closer to 1 cause the coefficients to change more quickly. Further, if it is assumed that within the relatively short cali- 
bration time, the noise source changes very little, i.e. SO n is equal to SO n _., in the equation above, these terms may 
40 be removed from the equation. 

[01 17] The second embodiment of the noise reduction system of the present invention will now be described in more 
detail. A high level block diagram of a second embodiment of the noise reduction system of the present invention 
utilizing a computational method to reduce echoes and oscillations is shown in Figure 11 . The noise reduction system 
of Figure 11, generally referenced 170, is constructed similarly to the noise reduction system 80 of Figure 4. The 
45 difference being the system 1 70 does not include the echo canceler circuit 95. The system 170 comprises an input 
transducer 1 72 such as a microphone, anti aliasing filter 174, fixed gain amplifier 176, A/D converter 178 and digital 
low pass filter 180. 

[01 1 8] The output of the low pass filter is corrected for non linearities of the transducers via a non linearity correction 
circuit comprising multiplier 1 86, sigma generator 1 83, coefficient processor 1 82 and linearity look up table (LUT) 1 84. 

50 [0119] The output of the multiplier is input to an input decoder 188 which functions to remove feedback picked up 
by the input transducer that was output by the output actuator. The output of the input decoder 188 is input to an 
equalizer 189 which comprises a digital filter that corrects the frequency response gain and group delay of the system 
analog elements including the output actuator and the input transducer. The result is that the frequency response of 
the combination of the input transducer and output actuator is flattened. The equalizer 189 causes the input signal, 

55 which lies within a particular frequency range, to be generated at the output transducer after a fixed delay. 

[0120] The output of the input decoder is amplified via a variable gain amplifier 192 whose gain is set by a gain 
control circuit 1 90. The output of the variable gain amplifier is input to the delayed cancellation circuit 1 94 whose output 
is added to the amplified signal via summer 1 96. The output of the summer is input to a D/A converter 1 97. The output 
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of the D/A converter is input to a low pass filter (LPF) 198. The output of the LPF is input to a power amplifier 200 
whose output drives the output actuator 202 which may comprise a low delay loudspeaker. 
[0121] Like components of the noise reduction system of Figure 11 function similarly to the corresponding compo- 
nents of the noise reduction system of Figure 4 and are thus described in more detail hereinabove. Note also that a 

5 hybrid I/O device can be used with the system 170 in place of a separate input transducer and output actuators. 

[0122] The system shown in figure 1 1 is based on the system time delay which is defined as the time from the output 
of the summer 186 through the input decoder 186, amplifier 192, summer 196, D/A 197, low pass filter 198, amplifier 
200, output actuator 202, acoustic medium, input transducer 1 72, filter 1 74, amplifier 1 76 and low pass filter 1 80. The 
time duration through the acoustic medium and the analog transducers and circuit comprises the time during which 

10 most of the energy of the impulse generated at the D/A converter 197 is measured at the A/D 1 78. The measurement 
(not shown) of the impulse is performed once during the start up phase of the system and considered as a constant 
during the system operation. 

[0123] The goal of the noise reduction system of Figure 11, as with the system of Figure 4 also, is to detect the 
relatively low noise signal from the larger signal picked up from the output of the loudspeakers. The performance of 

15 the system of Figure 11 is slightly less than the performance of the system of Figure 4 due to the fact that the echo 
canceler 95 functions to cancel acoustic echoes. The input decoder 1 88 functions to remove the effect of the feedback 
from the output actuators. Thus, the output of the decoder is substantially the noise signal with the output signal re- 
moved. This substantially pure noise signal is then inverted, equalized, amplified and output to the power amplifiers 
which drive the output actuators. As in the system of Figure 4, the delayed cancellation circuit functions to remove 

20 feedback that appears a system time delay later. Decoding of the input in the presence of a first or a second delayed 
cancellation signal that was fed back to the input requires that the delayed output be subtracted from each input signal 
sample. In addition, calculation of the linearity table in the presence of a first or a second delayed cancellation signal 
also requires that the delayed output be subtracted from each input signal sample. 

[0124] Another difference from the system of Figure 4 is that the digital data is output from the digital low pass filter 

25 1 80 to the input decoder 1 88 via the multiplier 1 86. The input decoder 1 88 functions to discern the interference noise 
signal from the input signal which includes a feedback signal having a relatively large amplitude. The input decoder 
comprises a £ generator similar to the £ generator 1 83 of the non linearity correction circuit and which is described in 
more detail below with reference to the flow diagram of Figure 13.. Note that the high level flow diagrams describing 
the method of linearity table calibration describe the calibration of a single coefficient. These method are repeated in 

30 order to calibrate all the coefficients. 

[0125] A high level flow diagram illustrating the first calibration method associated with the second embodiment is 
shown in Figure 12. The first method of calibration utilizes the fact that the noise is physical and continuous. The 
controller in the system tracks the relationship between values termed Table Input (Tl), Table Output (TO) and sigma 
(£) during operation of the system. The Tl values are measured at the input to the LUT 1 84. 

35 [01 26] The TO values are measured at the output of the multiplier 1 86. The £ values are generated by the £ generator 
183. The coefficient processor functions to calculate new LUT coefficients based on the Tl, TO and £ values. 
[0127] The calibration of the LUT coefficients during operation of the system attempts to ignore the effects of the 
noise source. Note that the input noise source itself changes between two adjacent samples.. The calibration method 
measures two adjacent TO values at the output of the multiplier 93 (Figure 11). Subsequently, the controller then 

40 measures two adjacent Tl values that are based on the previously measured TO values, assuming the system delay 
is known. The signal fed back from the output actuator to the input transducer is removed from the Tl value to yield an 
input value that reflects the noise level only. Subsequently, the effect of the noise source is removed from the Tl value. 
[0128] To aid in understanding the calibration method described herein, the method of removing feedback signals 
from the input signal implemented by the I generator 1 83 will first be described in more detail. A high level flow diagram 

45 illustrating the echo removal method of the present invention utilized in the second embodiment of the noise reduction 
system is shown in Figure 1 3. It is assumed that the initial speaker output is equal to zero and that the initial £ is equal 
to zero. The first step is to sample the input value and set £ to be equal to the input value in order to drive the output 
actuator 202 (Figure 11 ) (step 270). Next, the system waits for the effect of the output value to appear at the output of 
the LPF 180 , i.e., the Table Input (Tl) value, (step 272). The time needed for the output to propagate round trip is 

50 termed the system delay time. The system delay time must be known and is typically measured at the time the system 
is installed after it is first powered on. 

[0129] Once the effect of the output value appears at the Tl, the input value is read and denoted Tl r>1 (274). This 
input value, however, reflects both the noise level and the output that was fed back through the input. The system then 
starts examining input samples in order to detect a change in the input value. When a change is detected, a second 
55 input value is read and denoted Tl n (step 276). The delta A between the two input values is then calculated and added 
to a running sum £ (step 278). The £ value is then amplified in order to drive the actuator (step 280). The sum is 
initialized to zero at system startup time and is updated for each change in the input value. The running sum can be 
expressed as the following. 
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for the digital FIR filter that can be utilized in the delayed cancellation circuit. The system delay time is defined as the 
time from impulse generation until the time that the amplitude of the first response at the output of the A/D converter 
92 is equal to 0.2 of the maximum amplitude response. For example, if it is assumed that the system delay time is 1 00 
\is, then the first 100 us of the measured impulse response is discarded. The portion of the response 100 \ls and 
5 beyond is sampled and used to generate the filter coefficients. If the system includes cancellation circuitry for the third 
mirrored noise source, then the calibration is performed in steps. In the first step, the second mirrored source cancel- 
lation coefficient calibration is performed as described above. With this cancellation activated, calibration of the second 
step is performed by repeating the procedure described in Figure 5. The new coefficients are derived from the second 
feedback filter (not shown). 

10 [0071] The output of the digital summer 1 12 is input to the D/A converter 11 4 and to the echo canceler 95 portion of 
the system. The echo canceler 95 comprises a digital filter 96, an adaptation circuit 98 and a summer 94. The echo 
canceler 95 functions to remove reflections of sound and prevent oscillations caused by the output of the output actuato r 
feeding back to the input transducer. The filter 96 is a digital FIR filter having a sufficient number of taps, i.e., delay, to 
cover the round trip delay through the system as well as feed back echoes from the acoustical environment. 

is [0072] In order to reduce the computation burden of the FIR, interpolation and decimation methods, well known in 
the digital filter arts, are used (not shown) at the FIR output and input, respectively. The coefficients of the FIR filter 
correspond to the coefficients of the impulse response of the portion of the system from the output of the summer 1 1 2 
to the input of the summer 94 and are measured and setup at the time the system is started up. The measurements 
is performed by using the maximum length sequence technique (MLS), a technique well known in the art. The MLS 

20 technique is implemented by generating the MLS pseudo random sequence at the output of the summer 1 1 2, measuring 
at the input to summer 94, and then low pass filtering the cross correlation of the two signals. The results of the MLS 
measurement are the filter 96 coefficients. 

[0073] Real time adaptation of the digital filter coefficients by the adaptation circuit 98 is performed using any suitable 
FIR adaptation algorithm well known in the art, such as the least mean squares technique. The adaptation is based in 

25 the error signal output by the summer 94. 

[0074] Since the input noise signal from the noise source 24, e.g., motor or airflow noise, is a highly correlated signal, 
the regular FIR adaptation is applicable only outside of the auto correlation area. For example, in a typical air condi- 
tioning system, the auto correlation signal is 300 microseconds. Assuming a system delay of 1 00 microseconds, there 
are 300 minus 100 or 200 microseconds that may not be calibrated using the least mean square technique. 

30 [0075] If the length of the Fl R is adapted to be 5 milliseconds, then only the last 4.8 milliseconds may be adapted in 
real time using the least mean square (LMS) technique using a special variance explained bellow. Note that the LMS 
technique is well known in the art. The coefficients within the first 200 microseconds change relatively slowly with time 
and are effected mostly by temperature and humidity. Calibration of the coefficients within the first 200 microseconds 
is performed slowly in real time using the MLS technique described above with the exception that no generation of 

35 signal occurs. 

[0076] Instead, the noise data is used to perform cross correlation. Averaging a large number of test results yields 
the coefficients for the first 200 microseconds of the FIR filter 96. Note that when using the local closed loop control 
system to control the linearity of the output actuator as described in the linearity table, there is no need to adapt the 
first 200 microseconds in real time since the system response during that time is constant. 

40 [0077] The noise is cyclic, i.e., repetitive, and to enable acceptable adaptation, only one noise cycle is enabled within 
the length of the FIR filter. The full length of the FIR filer is divided into a plurality of short FIR filters. For example, if 
the noise cycle is 600 Hz, then the maximum allowed FIR length is 1/600 = 1.6 milliseconds. To cover the full 4.8 
milliseconds, three short FIR filters are used wherein each filter spans 1 .6 milliseconds. Each FIR filter comprises it's 
own separate adaptation circuit. 

45 [0078] The following example is regarded to LMS adaptation algorithm, but might be implemented on any known 
adaptation algorithm. The block diagram shown in Figure 6 illustrates the difference between a conventional LMS 
adaptation and the LMS adaptation of the present invention. Three FIR filter units are used with three taps in each FIR 
filter. The data from the summer 1 12 is input to the shift register 320 that functions as a serial to parallel converter. The 
first three registers 1 , 2, and 3 of the shift register 320 are utilized by the first FIR 321 in accordance with the following. 

50 The input noise that comes from multiplier 93 is multiplied by the constant u.1 via multiplier 326. The resulting product 
is multiplied with the contents of the first three registers of shift register 320 via multiplier 325. The resultant product 
is subsequently added to the current coefficient 323 via adder 324. 

[0079] It is important to note that the size of the coefficient 323 and the size of the FIR filter is 1/3 of the total length 
of the FIR filter. The remainder of the length is used as a shift register delay 322. The FIR portion 321 and the shift 
55 register delay portion 322 combine to form the complete FIR length. The summer 329 sums the output of the FIR with 
the input from multiplier 93 to generate first part of the result. The output of the summer is input to the summer 343. 
[0080] The second set of three registers 4, 5 and 6 are delayed in sync with delay 330 which generates a delay equal 
to the length of FIR 321 . The output of the delay 330 is input to FIR 331 . The coefficient 334 adaptation is performed 
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using multipliers 328 and 327 and adder 333 in a manner similar to that of the adaptation of coefficient 323 but utilizing 
constant u.2 and registers 4, 5 and 6. The output of the FIR 331 is input to delay 332 to complete the length of the FIR. 
The output of delay 332 is summed with the input from multiplier 93 via adder 335 to yield the second part of the result 
which is input to summer 343. 

5 [0081] The third portion of the output of summer 343 is generated utilizing bits 7, 8 and 9 of shift register 320 after 
passing the data through delay 340 and FIR 341 . The coefficient 339 adaptation is calculated using u.3, the input signal 
from multiplier 93, the output of multipliers 337, 336 and the output of adder 338. Note that the length of delay 340 is 
equal to the length of FIR 321 plus the length of FIR 331 . Thus, FIR 341 functions to complete the total length of the 
FIR filter. 

10 [0082] The dynamically changing values u.1 , u2 and u.3 function to determine how fast the adaptation algorithm 
performs. For example, if the adaptation algorithm runs too fast, the LMS will not yield correct results, i.e., it may not 
converge. A typical value for u.1 , u.2 and u.3 is in the range of 0.1 to 0.2 which was derived from experimentation. 
[0083] The output of the summer 94 is a digital signal wherein the feedback, i.e., echo, caused by the generated 
output signal has been removed. The digital filter 96 functions to generate a signal that is substantially the opposite of 

15 the portion of the input signal that is fed back from the output signal. 

[0084] The response of the adaptive digital filter 96 is the same response of that of the system from the output 
actuator to the input transducer including the environmental acoustic effects. The digital filter 96 functions to generate 
the signal output by the output transducer but with opposite polarity, i.e., negative polarity and the correct phase shift 
and amplitude. The summer 94 adds the output of the digital filter 96 with the output of the multiplier 93, which is the 

20 digital input corrected for the non linearities of the input transducer, to eliminate the echo from the input signal. The 
addition operation performed by summer 94 may be performed digitally as described above or can be performed in 
analog. The summing can be performed in analog by utilizing a D/A converter to convert the output of the digital filter 
96 to analog and summing the result with the analog signal from the input transducer. In the case of the analog echo 
cancellation summation, it is not economical to use an analog multiplier for the linearity table and thus, the assumption 

25 js that the actuator is linear. In either case, the output of the echo canceler is input to the digital low pass filter 1 00 as 
described hereinabove. 

[0085] It is important to note that the time response of the digital filter 96 must be long enough to cancel all of the 
acoustical sound fed back to the input transducer from the output actuator. The adaptive echo canceler 95 adjusts the 
filter taps in response to changes in the acoustical environment, e.g., temperature, furniture placement, etc. 

30 [0086] The output of the summer 112 is input to the D/A converter 114 which functions to convert the digital signal 
into analog. The output of the D/A converter 114 is input to an analog low pass filter 116. The cutoff frequency of the 
low pass filter 11 6 is calculated according to the D/A sampling rate and is designed to be low latency. The output of 
the low pass filter 1 1 6 is input to the power amplifier 1 1 8. The power amplifier boosts the level of the output sufficiently 
to drive the output actuator 120. The output actuator can be any suitable device that transforms the output signal to 

35 sound waves. The output actuator can comprise a loudspeaker having a relatively short delay. The output actuator 
120 may comprise a single output actuator or a plurality of output actuators connected in parallel. 
[0087] In addition , the system of Figure 4 also comprises a controller (not shown) which functions to administer and 
control the configuration, operation, settings and all timing of the noise reduction system 10. 
[0088] It is noted that the sound emitted from the output actuator is detected by the input transducer at a much higher 

40 amplitude than the noise itself. The noise reduction system 1 0 functions to discern the much smaller noise signal from 
the much larger echo signal. 

[0089] However, the input transducer exhibits a non linear response to the amplitude of the sound waves incident 
on it. It turns out that the input/output response of the input transducer, such as a microphone comprising a membrane, 
is dependent on the total amount of power hitting the microphone. Thus, for the same amplitude noise signal, the 

45 linearity of the microphone is different for an input with low total power input as compared to an input with high total 
power. Thus, the high amplitude signal output from the output actuator effects a different relative output signal from 
the input transducer due to the noise. This is especially true when using a hybrid to perform both the input and the 
output transducer functions of the system. Non linearity at the output actuator is manifested when the output peak 
power to the transducer is large. When using separate input/output devices, the non linearity in the output device is 

50 much more critical when the mechanical amplitude is large. The operation of the system is based on the correct dis- 
crimination of the noise source from the much larger echo signal. Thus, a correction of the echo signal amplitude as 
it appears in the input is required in order to achieve good cancellation of the echo signal. The linearity LUT performs 
this correction operation by multiplying each input value with the correction coefficient. 

[0090] The linearity table may be generated a priori and calibrated at the time of installation or adjusted dynamically 
55 on the fly. The linearity table and the dynamic adjustment of its contents are described in more detail hereinbelow. 
[0091] As an alternative to the linearity table, a local closed loop control system can be used to compensate for the 
linearity of the actuator (not shown). Such a technique is described, for example, in the proceedings of Active 97 in 
the article by Yoon-Sun Kim and Youngjin Park entitled "Non Linearity Compensation for Harmonic Distortion of Direct 
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Radiation Loudspeaker." 

[0092] A high level block diagram illustrating the input transducer and output actuator implemented as a hybrid com- 
bination in a single element is shown in Figure 7. As an alternative to a separate input transducer and output actuator, 
a hybrid input/output (I/O) element or transducer 1 62 can be utilized that performs the functions of both elements. 

5 [0093] In this case, a hybrid circuit 1 60 is needed to interface the I/O element 1 62 to the power amplifier 1 64 and 
the anti aliasing filter 1 66. The hybrid 1 60 functions to transfer power from the power amplifier 1 64 with minimum losses 
to the I/O transducer element 162 and with maximum reduction to the anti aliasing filter 1 66. 
[0094] The gain control method implemented in the gain control circuit 1 06 (Figure 4) of the variable gain amplifier 
108 will now be described in more detail. A high level flow diagram illustrating the gain control method utilized in both 

w the first and second embodiments is shown in Figure 8. This method operates in a loop to continuously search for 
oscillations in the system. In general, an FFT or DFT performed in the gain control circuit is used to map the mean 
amplitude of the frequency content in the detected input. When a new frequency is detected or when the total input 
power to the system increases, an immediate reduction of the gain of the amplifier 108 is performed. 
[0095] The first step is to check the total power of the output signal to see if it is over a predetermined maximum 

15 (step 1 30). If the power is over the permitted maximum then the gain is reduced until the total output power is below 
the maximum (step 132). Once the gain is within the permitted range, the current gain setting is stored as the initial 
gain value (step 134). A Fast Fourrier Transform (FFT) or DFT is then performed on the signal input to the amplifier 
(step 1 36). The FFT or DFT yields a map of the frequency content of the input signal. A plurality of samples of the input 
are taken and FFT or DFT analysis performed on each sample. Corresponding frequency elements are averaged over 

20 many samples (step 138). 

[0096] It is then checked for the significant presence of signal content at new frequencies (step 1 40). If signal content 
is found at new frequencies, these frequencies elements are tracked over time (step 142) to determine whether they 
are oscillations (step 144). A new frequency that is persistent in time is considered a suspect oscillation. If there is a 
suspected oscillation, then the gain is reduced by a step amount (step 146). If the gain is reduced to a predetermined 

25 minimum (step 148) the frequencies that were suspect as oscillations are mapped as input noise, i.e., a frequency 
generated by the environment with a particular amplitude and frequency range and not as oscillations (step 152). This 
is because, at such low values of gain, it would be highly improbable that the suspect frequencies were caused by 
oscillations. Signals within the frequency range are treated as environmental noise. After the environmental frequencies 
are mapped, the gain is restored to the gain that was initially saved during step 134. 

30 [0097] If there are no new frequency elements in step 140 or oscillations suspected in step 144, the method ends. 
If the gain in step 148 is not a minimum, than another group of samples is input (step 150). After sampling another 
group of input signals, it is then checked whether oscillations are still present after the gain has been reduced. The 
method keeps looping by checking for oscillations, reducing the gain if oscillations are found and generating additional 
input sample until either no oscillations are found or the gain reaches a minimum. 

35 [0098] As described previously, the linearity LUT 97 is used to correct the non linear behavior of the input/output 
transducers. Since typical input transducers are constructed from moving mechanical parts, the output response for 
the same input noise level is different for different levels of total input power incident on the input transducer. However 
the output actuator suffers from non linearities when the motion of the moving mechanical parts increase as a result 
of increasing output power. The difference in output response of the output actuator and input transducer is corrected 

40 utilizing the linearity LUT 97. 

[0099] The linearity LUT stored different coefficients for each amplitude input value range. The linearity LUT is divided 
into regions with each region having its own coefficient value. For example, assuming a 12 bit data word, the linearity 
LUT may be divided into 256 regions, thus only utilizing 8 bits of the input data. Upon startup of the noise reduction 
system 10, all of the coefficients in the linearity LUT have the same value. 

45 [01 00] The present invention includes a first and a second method of calibrating the linearity LUT. Both methods are 
performed in real time during the operation of the noise reduction system. Both calibration methods utilize a linearity 
look up table (LUT) that holds coefficients used to adjust the input data output by the A/D converter. The adaptation 
or calibration of the coefficients of the linearity LUT serves to compensate for slow changes to the linearity of the system 
caused by temperature, humidity etc. The linearity LUT also functions to compensate forthe mechanical non linearities 

50 of both the input transducer and the output actuator(s). Although not necessary to perform the invention, the input 
amplitude is divided into a plurality of regions wherein each region has a coefficient associated with it. 
[0101] The first calibration method utilizes the input noise signal itself to calibrate the coefficients. In the second 
method, a small amplitude calibration signal is injected into the system du ring operation and the results used to generate 
new coefficient. Each calibration method will now be described in more detail below. The high level flow diagrams 

55 describing the linearity table calibration methods refer to the calibration of one coefficient. The methods described are 
repeated in order to calibrate all the coefficients. 

[01 02] A high level flow diagram illustrating the first calibration method associated with the first embodiment is shown 
in Figure 9. The first method of calibration utilizes the fact that the noise is physical and continuous. The controller in 
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the system tracks the relationship between values termed Table Input (Tl), Table Output (TO) and Summer Output 
(SO) during operation of the system. The Tl values are measured at the output of the A/D converter 92 (Figure 4), the 
TO values are measured at the output of the multiplier 93 and the SO values are measured at the output of the summer 
94. The coefficient processor 99 functions to calculate new LUT coefficients based on the Tl, TO and SO values. 
5 [0103] The calibration of the LUT coefficients during operation of the system attempts to ignore the effects of the 
noise source. Note that the input noise source itself changes between two adjacent samples. Note also that the output 
of the summer SO represents the noise source since the echo canceler 95 effectively removes the echo signal from 
the input signal. 

[0104] The following calibration method is based on the system time delay which is defined as the time from TO to 
10 the D/A 114 plus the time from the D/A 114 to the A/D 92. Note that factored in this time is the acoustic medium, the 
analog transducers the circuit and the time from the A/D 92 to TO. The time attributed to the acoustic medium and the 
analog transducers and circuit is substantially equivalent to the time during which most of the energy of the impulse 
generated at D/A 11 4 is measured at the A/D converter 92. 

[0105] The measurement of the impulse is performed using the impulse response measurement technique described 

15 previously in connection with the calibration of FIR filter 96. 

[0106] The calibration method initially measures two TO values, i.e., a value at time 'n-1' and another at time 'n\ 
Subsequently, after a system time delay, the controller then measures the two Tl and SO values that are the effect of 
the previously measured TO values, assuming the system time delay is known. In other words, each Tl and SO value 
is measured after the corresponding TO values has had a chance to propagate through the system, i.e., output by the 

20 output actuator, sensed at the input transducer, etc. Note that the SO value is the Tl value after compensation for non 
linearities and after the echo fed back from the output actuator to the input transducer is removed to yield a value that 
reflects the noise level only. 

[0107] With reference to Figure 9, the first step in the calibration method is to read a Table Output (TO) value at the 
output of the multiplier 93 (Figure 4), termed TO n . 1 (step 21 0). After reading the TO^ value, the system waits for the 

25 output of the multiplier to change before reading the next value which is termed TO n (step 21 4). During operation of 
the system, it is not unusual if the input data does not vary much from sample to sample. This is due to the fact that 
the sampling rate of the system is orders of magnitude higher than the frequencies making up the noise signal. In step 
214, the system waits more then one sample time, and the sample is taken when the difference between the level of 
TO^ and TO n reaches some minimum which is predetermined. This minimum determines the coefficient calibration 

30 accuracy. It is then determined whether TO n is within the same region of the LUT as TO,,.., (step 21 6). This is checked 
in order to prevent two TO values being associated with different regions of the LUT. The calibration method calculates 
new coefficients for a single region of the LUT at a time. The calculations, thus, cannot span borders between regions. 
[0108] Next, the system waits for the effect of the TO^ value (step 21 8) and TO n value (step 222) to appear at the 
output of the A/D converter an d the outp ut of the summer . As described previously, the data output by the A/D converter 

35 is termed Table Input (Tl) data and the data output by the summer is termed Summer Output (SO). Once the effect of 
the TO^ data appears at the output of the A/D converter and the summer, the Tl rv1 and SO n _., values are read (step 
220). Similarly, once the effect of the TO n data appears at the output of the A/D converter and the summer, the Tl n and 
the SO n values are read N n (step 224). The steps of first waiting and then reading the Tl and SO values described 
herein can be implemented either sequentially or in parallel. 

40 [0109] The index to the LUT which determines which coefficient is presently under calculation is then generated 
using the Tl n .., value read during step 220 (step 225). This index is used for the calibration process only and does not 
effect the main real time data path of the system. In the example noise reduction system shown in Figure 4, the LUT 
has less entries in it than the number of possible input values, e.g., 256 regions for 12 bits of input data, in order to 
reduce the size of the lookup table required. Alternatively, the LUT can be constructed to hold a coefficient value for 

45 each and every possible input data. 

[0110] Once the Tl and SO data have been read, the new coefficients are calculated using the following equation 
(step 226). 



50 



Cnew = Cold + K ] 



TO -TO.. 



(27. -SO.) -(27^- SO..^ 



55 [0111] The new coefficient C new is a function of the old coefficient C old . The values TO, Tl and SO are used to 
generate an intermediate new coefficient from which C ojd is subtracted. A portion of the delta is added to C oW to perform 
the calibration. The constant K varies between 0 and 1 and is used to determine the speed with which the coefficients 
are permitted to change. Values of K closer to 0 cause the coefficients to change more slowly whereas values of K 
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closer to 1 cause the coefficients to change more quickly. 

[0112] The second calibration method associated with the first noise reduction embodiment of Figure 4 will now be 
described in more detail. A high level flow diagram illustrating the second calibration method associated with the first 
embodiment is shown in Figure 10. This second method of calibration is similar to that of the first method described in 

5 connection with Figure 9, with the difference being that rather than wait for the actual noise source to cause a change 
to the TO value, an artificial noise signal is injected into the data path to simulate a known change in the noise signal 
level. The second method of coefficient calibration models the noise as a pseudo DC level noise source during cali- 
bration. This is a reasonable assumption since the calibration period of approximately 10 \is is very short relative to 
the noise frequency. This pseudo DC level of the noise is used to point to a particular region in the linearity LUT 

10 [01 1 3] With reference to Figures 4 and 1 0, the first step is to measure the Table Input (Tl) value denoted Tl^ at the 
output of the A/D converter 92 (step 230). The index to the linearity LUT is then generated based on the Tl^ value 
just measured (step 232). The index determines which of the coefficients of the linearity LUT is to be calibrated during 
this particular invocation of the method. The next step is to read the Table Output (TO) value, denoted TO^, at the 
output of the multiplier 93 and the Summer Output (SO) value, denoted SO,,..,, at the output of the summer 94 (step 

15 234). The TO value is generated by multiplying the output of the A/D converter with the output of the LUT The result 
of the multiplication is added to the output of the digital filter. Note that the TO and SO values are read immediately 
after the Tl value is read without waiting a system time delay as in the first method of Figure 9. 
[0114] A calibration signal is then injected at the output of the multiplier (step 236). The output of the multiplier which 
is denoted as the TO value is replaced with the calibration signal for a finite time period. The calibration signal, termed 

20 TO n , comprises the original output of the multiplier TO n _., increased by a known delta amount. The system then waits 
one system delay time for the injected calibration signal to appear at the output of the A/D converter (step 238). After 
waiting one system time delay, the data at the output of the A/D converter is read and termed Tl n . In addition, the output 
of the summer 94, termed SO n , is also read (step 240). 

[01 1 5] After reading the Tl n value, the original output of the multiplier TO n _., before the calibration signal was injected 
25 js restored (step 242). Then, based on the values Tl, TO and SO, the new coefficient can be calculated utilizing the 
following equation (step 244). 



30 



C nnu = C old + K 



TO, - TO. 



p tt -SoJi-(TC-SO n _ { ) 



-c, 



[0116] The new coefficient C new is a function of the old coefficient C okj . The values TO, Tl and SO are used to 
35 generate an intermediate new coefficient from which C 0 , d is subtracted. A portion of the delta is added to C 0 , d to perform 
the calibration. The constant K varies between 0 and 1 and is used to determine the speed with which the coefficients 
are permitted to change. Values of K closer to 0 cause the coefficients to change more slowly whereas values of K 
closer to 1 cause the coefficients to change more quickly. Further, if it is assumed that within the relatively short cali- 
bration time, the noise source changes very little, i.e. SO n is equal to SO n .., in the equation above, these terms may 
40 be removed from the equation. 

[0117] The second embodiment of the noise reduction system of the present invention will now be described in more 
detail. A high level block diagram of a second embodiment of the noise reduction system of the present invention 
utilizing a computational method to reduce echoes and oscillations is shown in Figure 11 . The noise reduction system 
of Figure 11, generally referenced 170, is constructed similarly to the noise reduction system 80 of Figure 4. The 
45 difference being the system 170 does not include the echo canceler circuit 95. The system 170 comprises an input 
transducer 172 such as a microphone, anti aliasing filter 174, fixed gain amplifier 176, A/D converter 178 and digital 
low pass filter 180. 

[01 1 8] The output of the low pass filter is corrected for non linearities of the transducers via a non linearity correction 
circuit comprising multiplier 1 86, sigma generator 1 83, coefficient processor 1 82 and linearity look up table (LUT) 1 84. 

so [0119] The output of the multiplier is input to an input decoder 188 which functions to remove feedback picked up 
by the input transducer that was output by the output actuator. The output of the input decoder 188 is input to an 
equalizer 1 89 which comprises a digital filter that corrects the frequency response gain and group delay of the system 
analog elements including the output actuator and the input transducer. The result is that the frequency response of 
the combination of the input transducer and output actuator is flattened. The equalizer 189 causes the input signal, 

55 which lies within a particular frequency range, to be generated at the output transducer after a fixed delay. 

[0120] The output of the input decoder is amplified via a variable gain amplifier 192 whose gain is set by a gain 
control circuit 1 90. The output of the variable gain amplifier is input to the delayed cancellation circuit 1 94 whose output 
is added to the amplified signal via summer 1 96. The output of the summer is input to a D/A converter 1 97. The output 
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of the D/A converter is input to a low pass filter (LPF) 198. The output of the LPF is input to a power amplifier 200 
whose output drives the output actuator 202 which may comprise a low delay loudspeaker. 
[0121] Like components of the noise reduction system of Figure 11 function similarly to the corresponding compo- 
nents of the noise reduction system of Figure 4 and are thus described in more detail hereinabove. Note also that a 

5 hybrid I/O device can be used with the system 1 70 in place of a separate input transducer and output actuators. 

[0122] The system shown in figure 1 1 is based on the system time delay which is defined as the time from the output 
of the summer 1 86 through the input decoder 1 86, amplifier 1 92, summer 1 96, D/A 1 97, low pass filter 1 98, amplifier 
200, output actuator 202, acoustic medium, input transducer 1 72, filter 1 74, amplifier 1 76 and low pass filter 1 80. The 
time duration through the acoustic medium and the analog transducers and circuit comprises the time during which 

10 most of the energy of the impulse generated at the D/A converter 197 is measured at the A/D 178. The measurement 
(not shown) of the impulse is performed once during the start up phase of the system and considered as a constant 
during the system operation. 

[0123] The goal of the noise reduction system of Figure 11 , as with the system of Figure 4 also, is to detect the 
relatively low noise signal from the larger signal picked up from the output of the loudspeakers. The performance of 

15 the system of Figure 11 is slightly less than the performance of the system of Figure 4 due to the fact that the echo 
canceler 95 functions to cancel acoustic echoes. The input decoder 1 88 functions to remove the effect of the feedback 
from the output actuators. Thus, the output of the decoder is substantially the noise signal with the output signal re- 
moved. This substantially pure noise signal is then inverted, equalized, amplified and output to the power amplifiers 
which drive the output actuators. As in the system of Figure 4, the delayed cancellation circuit functions to remove 

20 feedback that appears a system time delay later. Decoding of the input in the presence of a first or a second delayed 
cancellation signal that was fed back to the input requires that the delayed output be subtracted from each input signal 
sample. In addition, calculation of the linearity table in the presence of a first or a second delayed cancellation signal 
also requires that the delayed output be subtracted from each input signal sample. 

[0124] Another difference from the system of Figure 4 is that the digital data is output from the digital low pass filter 

25 1 80 to the input decoder 1 88 via the multiplier 1 86. The input decoder 1 88 functions to discern the interference noise 
signal from the input signal which includes a feedback signal having a relatively large amplitude. The input decoder 
comprises a Z generator similar to the I generator 1 83 of the non linearity correction circuit and which is described in 
more detail below with reference to the flow diagram of Figure 13.. Note that the high level flow diagrams describing 
the method of linearity table calibration describe the calibration of a single coefficient. These method are repeated in 

30 order to calibrate all the coefficients. 

[0125] A high level flow diagram illustrating the first calibration method associated with the second embodiment is 
shown in Figure 12. The first method of calibration utilizes the fact that the noise is physical and continuous. The 
controller in the system tracks the relationship between values termed Table Input (Tl), Table Output (TO) and sigma 
(Z) during operation of the system. The Tl values are measured at the input to the LUT 184. 

35 [01 26] The TO values are measured at the output of the multiplier 1 86. The Z values are generated by the Z generator 
183. The coefficient processor functions to calculate new LUT coefficients based on the Tl, TO and Z values. 
[0127] The calibration of the LUT coefficients during operation of the system attempts to ignore the effects of the 
noise source. Note that the input noise source itself changes between two adjacent samples.. The calibration method 
measures two adjacent TO values at the output of the multiplier 93 (Figure 11). Subsequently, the controller then 

40 measures two adjacent Tl values that are based on the previously measured TO values, assuming the system delay 
is known. The signal fed back from the output actuator to the input transducer is removed from the Tl value to yield an 
input value that reflects the noise level only. Subsequently, the effect of the noise source is removed from the Tl value. 
[0128] To aid in understanding the calibration method described herein, the method of removing feedback signals 
from the input signal implemented by the Z generator 1 83 will first be described in more detail. A high level flow diagram 

45 illustrating the echo removal method of the present invention utilized in the second embodiment of the noise reduction 
system is shown in Figure 1 3. It is assumed that the initial speaker output is equal to zero and that the initial Z is equal 
to zero. The first step is to sample the input value and set Z to be equal to the input value in order to drive the output 
actuator 202 (Figure 11) (step 270). Next, the system waits for the effect of the output value to appear at the output of 
the LPF 180 , i.e., the Table Input (Tl) value, (step 272). The time needed for the output to propagate round trip is 

50 termed the system delay time. The system delay time must be known and is typically measured at the time the system 
is installed after it is first powered on. 

[0129] Once the effect of the output value appears at the Tl, the input value is read and denoted Tl^ (274). This 
input value, however, reflects both the noise level and the output that was fed back through the input. The system then 
starts examining input samples in order to detect a change in the input value. When a change is detected, a second 
55 input value is read and denoted Tl n (step 276). The delta A between the two input values is then calculated and added 
to a running sum Z (step 278). The Z value is then amplified in order to drive the actuator (step 280). The sum is 
initialized to zero at system startup time and is updated for each change in the input value. The running sum can be 
expressed as the following. 



14 



EP 1 414 021 A1 



I = Z + A 

[0130] These calculations are performed each time there is a change in Tl except during the period within a system 
delay and one clock period after system delay, as shown in the Table 1 below. As an alternative, the above calculations 
may be performed once every system clock except during the period that is during the system delay and one clock 
10 after system delay, as shown in Table 1 below, regardless of the resulting value of A. Note that the value of A may often 
be zero due to the slowly varying input signal relative to the sampling frequency of the A/D converter. 
[01 31 ] The following table is presented to illustrate the method of Figure 1 3. 
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[0132] The following notes will aid in understanding the contents of Table 1 presented above. 

1 . Note that in system operation period 1 Tl^ has no value, thus according to the equation given above, A and Z 
are not defined. In this case, Z is set equal to the input value. 
5 2. The 'X' in the A calculations row means that no A calculations are performed in that particular system operation 

period. 

3. The time difference between system operation periods that are not marked with '10 \is delay' is time between 
samples taken by the A/D converter 178 (Figure 11). Note that there may be numerous periods between 'system 
delay' periods, i.e., numerous periods until the noise input value changes. 
10 4. Since the system functions to accumulate errors in measurement accuracy, the Z accumulator does not accu- 

mulate continuously. Periodically, before the error increases above a predetermined value, the system waits until 
Z becomes nearly zero at which time it sets the Z generator to the initial conditions. Thus, regardless of the history, 
the system starts over from period 1 as shown in Table 1 . 

5. The contents of Table 1 above are associated with the operations of the I generator 1 83 and the input decoder 
15 1 88 (Figure 1 1 ). It is preferable to set the operation of the two generators to the initial state as explained in Note 

4 at the substantially the same time. 

[0133] The total measured input is the sum of the noise input and the speaker output. The A and Z values are 
calculated as given above and the speaker output represents the value output by the output actuator as sensed in the 
20 input. In this example, the system is considered to have a gain of 1 0, thus the speaker output, as fed back to the input, 
is taken as 1 0 times the value of Z. As can be seen from Table 1 above, the method effectively removes the effect of 
the speaker output from the input data. 

[0134] With reference to Figures 11 and 12, the first step in the calibration method is to read a Table Output (TO) 
value at the output of the multiplier 182, termed TO^ (step 300). 
25 [0135] After reading the TO n .-, value, the system waits for the difference between TO n _., and TO n at the output of the 
multiplier to reach a predetermined minimum before reading the next TO value which is termed TO n (step 302). This 
minimum determines the accuracy of the coefficient calibration. 

[0136] It is then determined whetherTO n is within the same region of the LUT as TO n _., (step 304). This is checked 
in order to prevent two TO values being associated with different regions of the LUT The calibration method calculates 

30 new coefficients for a single region of the LUT at a time. The calculations, thus, cannot span borders between regions. 
[0137] Next, the system waits for the effect of the TO n .., value (step 306) and TO n value (step 31 0) to appear at the 
output of the LPF. The data output by the LPF is termed Table Input (Tl) data. Once the effect of the TO rv1 data appears 
at the output of the LPF, the Tl n .., value is read along with the value of the noise denoted N n .., (step 308). N is the net 
noise extracted from the echo removal method described hereinabove and is equal to the Z at that particular point in 

35 time. Similarly, once the effect of the TO n data appears at the output of the LPF, the Tl n value is read along with the 
value of the noise denoted N n (step 312). The steps of first waiting and then reading the Tl values described above 
can be implemented either sequentially or in parallel. 

[0138] The Tl n . 1 value is then used to generate an index into the LUT (step 314). The LUT has less entries in it than 
the number of possible input values, e.g., 256 regions for 12 bits of input data in order to reduce the size of the LUT. 
40 Alternatively, the LUT can be constructed to hold a coefficient value for each and every possible input data. 

[0139] Once the Tl and N data have been read, the new coefficient is calculated using the following equation. 



45 



TO-TO K _ x 



-c, 



old I 



[01 40] The new coefficient C new is a function of the old coefficient Cold. The values TO, Tl and Z are used to generate 
50 an intermediate new coefficient from which the old coefficient C old is subtracted. A portion of the delta (determined by 
the constant K) is added to C old to perform the calibration. The constant K varies between 0 and 1 and is used to 
determine the speed with which the coefficients are permitted to change. Values of K closer to 0 cause the coefficients 
to change more slowly whereas values of K closer to 1 cause the coefficients to change more quickly. 
[0141] The second calibration method associated with the first noise reduction embodiment of Figure 4 will now be 
55 described in more detail. A high level flow diagram illustrating the second calibration method associated with the second 
embodiment of the present invention is shown in Figure 14. This second method of calibration is similar to that of the 
first method described in connection with Figure 12, with the difference being that rather than wait for the actual noise 
source to cause a change to the TO value, an artificial noise signal is injected into the data path to simulate a known 
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change in the noise signal level. The method utilizes the output of the LPF 1 80 (Tl values), the output of the multiplier 
1 86 (TO values) and the output of the I generator 1 83 (I values) in performing the calibration calculations. The second 
method of coefficient calibration is performed very quickly relative to the frequency of the input noise. In addition, it is 
assumed that the smallest period of the input noise is small enough relative to the time delay of the system that it can 

5 be regarded as a DC level. This pseudo DC level of the noise is used to point to a particular region in the linearity LUT. 
[0142] With reference to Figures 11 and 14, the first step is to measure the Table Input (Tl) value denoted Tl^ at 
the output of the LPF. The output of the Z generator is also calculated and denoted I n .-, (step 250). The index to the 
linearity LUT is then generated based on the Tl^ value just measured (step 252). The index determines which of the 
coefficients of the linearity LUT is to be calibrated during this particular invocation of the method. The next step is to 

10 read the Table Output (TO) value, denoted TO^ , at the output of the multiplier 1 86 (step 254). The TO value is gen- 
erated by multiplying the output of the LPF with the output of the LUT. The result of the multiplication is input to the 
input decoder. Note that the TO value is read immediately after the Tl value is read without waiting a system time delay. 
[0143] A calibration signal is then injected at the output of the multiplier (step 256). The output of the multiplier which 
is denoted as the TO value is replaced with the calibration signal for a finite time period. The calibration signal, termed 

15 TO n , comprises the original output of the multiplier TO n _i increased by a known delta amount. The system then waits 
one system delay time for the injected calibration signal to appear at the output of the LPF (step 258). After waiting 
one system time delay, the data at the output of the LPF is read and termed Tl n . In addition, the output of the I generator 
183 is read and termed £ n (step 260). 

[01 44] After reading the Tl n value, the original output of the multiplier TO n ..| before the calibration signal was injected 
20 is restored (step 262). Then, based on the values Tl, TO and I, the new coefficient is calculated utilizing the following 
equation (step 264). 
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The new coefficient C new is a function of the old coefficient C oJd . The values TO, Tl and Z are used to generate an 
30 intermediate new coefficient from which C o)d is subtracted. A portion of the delta is added to C okJ to perform the cali- 
bration. The constant K varies between 0 and 1 and is used to determine the speed with which the coefficients are 
permitted to change. Values of K closer to 0 cause the coefficients to change more slowly whereas values of K closer 
to 1 cause the coefficients to change more quickly. 

[0145] Note that since the calibration calculation for each coefficient occurs very quickly, the input noise source can 
35 be regarded as constant or a pseudo DC value at the time the Tl values are measured. Thus, the Z n _., and 1^ values 
will cancel and thus can be removed from the equation. 

[01 46] While the invention has been described with respect to a limited number of embodiments, it will be appreciated 
that many variations, modifications and other applications of the invention may be made. 

40 

Claims 

1 . An acoustic noise reduction system for reducing the effects of a noise source, comprising: 

45 input transducer means for sensing the acoustic noise field generated by the noise source and for generating 

an input signal therefrom; 

output actuator means for generating an acoustic output field that is effective to reduce the level of the acoustic 
noise field; 

correction means for adjusting the input signal generated by said input transducer to compensate for the non 
so linear characteristics of said input transducer and output actuator; 

echo cancellation means for removing from the input signal a portion of the output of said output actuator 
means fed back through said input transducer means, the output of said echo cancellation means representing 
a signal corresponding to substantially the noise source by itself; 

antinoise means for generating an antinoise signal opposite in phase to said input signal, said output actuator 
55 means generating said acoustic output field from said antinoise signal; and 

wherein said input transducer means is located in relatively close proximity to said output actuator means. 
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2. An acoustic noise reduction system for reducing the effects of a noise source, comprising: 

input transducer means for sensing the acoustic noise field generated by the noise source and for generating 
an input signal therefrom; 

output actuator means for generating an acoustic output field that is effective to reduce the level of the acoustic 
noise field; 

correction means for adjusting the input signal generated by said input transducer to compensate for the non 
linear characteristics of said input transducer; 

input decoding means for removing extraneous signals from said input signal so as to generate a signal cor- 
responding to substantially the noise source alone; 

antinoise means for generating an antinoise signal opposite in phase to said input signal, said output actuator 
means generating said acoustic output field from said antinoise signal; and 

wherein said input transducer means is located in relatively close proximity to said output actuator means. 

3. The system according to claims 1 , 2, 32 or 33 wherein said input transducer means comprises a microphone. 

4. The system according to claims 1 , 2, 32 or 33 wherein said output actuator means comprises at least one loud- 
speaker. 

5. The system according to claim 1 or 32, wherein said correction means comprises: 

storage means for storing a plurality of coefficients; 

coefficient processing means for dynamically updating the values of said plurality of 
coefficients stored in said storage means; and 

means for generating a corrected input signal from the contents of said storage means and said input signal. 

6. The system according to claim 2 or 33, wherein said correction means comprises: 

storage means for storing a plurality of coefficients; 

sigma generating means for outputting a signal corresponding to substantially the noise source only; 
coefficient processing means for dynamically updating the values of said plurality of coefficients stored in said 
storage means; and 

means for generating a corrected input signal from the contents of said storage means and said input signal. 

7. The system according to claims 5 or 6, wherein said means for generating a corrected input signal comprises a 
multiplier. 

8. The system according to claims 5 or 6, wherein said storage means is divided into regions whereby a particular 
range of input signal values map to the same coefficient value. 

9. The system according to claims 1, 2, 32 or 33, further comprising calibration means for dynamically updating said 
correction means in response to the input signal and the output of said correction means. 

10. The system according to claims 1 , 2, 32 or 33, further comprising an equalizer for compensating the frequency 
response gain and group delay of said system. 

11. The system according to claim 1 or 32, wherein said echo cancellation means comprises: 

a digital filter having a delay line with a number of taps whose total delay time is equivalent to at least a system 
time delay of said noise reduction system; 

adaptation means for dynamically adjusting the coefficient values associated with each of the taps of said 
digital filter; and 

summing means for adding the output of said digital filter with the output of said correction means. 

12. The system according to claim 1 or 32, wherein said echo cancellation means comprises: 

a shift register whose parallel outputs are divided into a plurality of N portions; 
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a plurality of N finite impulse response (FIR) means wherein each N m portion of said shift register coupled to 
the input of the N m FIR means; 

a plurality of N adaptation circuits wherein one adaptation circuit is associated with each of N FIR means; and 
wherein said N FIR means combine to yield an FIR filter whose length is equivalent to the combined length 
5 of said N FIR means, whereby the length of each FIR means is less then or equal to the period of a noise cycle. 

13. The system according to claims 1 , 2, 32 or 33, wherein said antinoise means comprises; 

a variable gain amplifier operative to generate an amplified signal 180 degrees opposite in phase from said 
10 input signal; and 

gain control means for dynamically controlling the gain of said variable gain amplifier. 

14. The system according to claim 13, wherein said gain control means is adapted to receive a manual input control 
signal from a user which determines the gain of said variable gain amplifier, said user able to vary the location of 

15 a quiet zone generated by said system by varying said input control signal. 

15. The system according to claim 14, wherein said input control signal is generated by said user remotely from said 
system and transmitted to said system via wireless communication means. 

20 1 6. The system according to claims 1 , 2, 32 or 33, further comprising anti alias filter means for removing alias frequen- 
cies from said input signal. 

17. The system according to claim 16, wherein said anti alias filter means comprises a low pass filter having a cutoff 
frequency sufficiently high enough so as to reject frequencies greater then the sampling rate utilized within said 

25 system. 

18. The system according to claims 1 , 2, 32 or 33, further comprising a low pass filter operative to reduce oscillations 
present in the system derived from feedback of the acoustic output field to said input transducer. 

30 19. The system according to claims 1 , 2, 32 or 33, further comprising delay cancellation means for reducing the effect 
of echo signals caused by said antinoise means sensed by said input transducer. 

20. The system according to claim 1 9, wherein said delay cancellation means comprises a digital filter whose output 
is added to the output of said antinoise means, to cancel speakers echo. 

35 

21 . The system according to claim 1 9, wherein said delay cancellation means comprises a plurality of delay cancellation 
circuits wherein each delay cancellation circuit is operative to reduce the effect of the echo caused by previous 
delay cancellation circuits. 

40 22. The system according to claim 20, wherein said digital filter comprises a finite impulse response (FIR) digital filter. 

23. The system according to claims 1 , 2, 32 or 33, further comprising means for sampling the acoustic noise field at 
a sampling rate approximately 1000 times or greater than the frequency of the noise source. 

45 24. The system according to claims 1 , 2, 32 or 33, wherein said system is operative to generate a controlled far field 
quiet zone. 

25. The system according to claims 1 , 2, 32 or 33, further comprising means for sampling a particular point of the 
noise source having an arbitrary phase and effectively canceling noise radiated from other points of the noise 

50 source having the same phase. 

26. The system according to claims 1 , 2, 32 or 33, wherein said input transducer means comprises a single microphone, 
said output actuator means comprises a plurality of loudspeakers arranged so as to radiate energy over a wide 
angle, the superimposed energy of said loudspeakers creates the total energy needed to control the noise over 

55 said wide angle. 

27. The system according to claims 1 , 2, 32 or 33, wherein said input transducer means comprises a single microphone, 
said output actuator means comprises a plurality of loudspeakers arranged in linear fashion, each loudspeaker 
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radiating energy in the same direction so as to create a wall of silence. 

28. The system according to claims 1, 2, 32 or 33, wherein said input transducer means and said output actuator 
means are implemented as a single input/output (I/O) hybrid device. 

5 

29. The system according to claims 1 , 2, 32 or 33, further comprising means for generating an audio source from said 
output actuator means. 

30. A method for reducing the effects of a noise source, comprising the steps of: 

10 

sensing the acoustic noise field generated by the noise source and generating an input signal therefrom; 
generating an acoustic output field that is effective to reduce the level of the acoustic noise field; 
adjusting the input signal generated by an input transducer to compensate for the non linear characteristics 
of said input transducer and an output actuator; 
15 removing from the input signal a portion of the output of the output actuator fed back through said input trans- 

ducer; 

generating a signal corresponding to substantially the noise source by itself; and 

generating an antinoise signal opposite in phase to said input signal, generating said acoustic output field from 
said antinoise signal. 

20 

31. A method for reducing the effects of a noise source, comprising the steps of: 

sensing the acoustic noise field generated by the noise source and for generating an input signal therefrom; 
generating an acoustic output field that is effective to reduce the level of the acoustic noise field; 
25 adjusting the input signal generated by an input transducer to compensate for the non linear characteristics 

of said input transducer; 

removing extraneous signals from said input signal so as to generate a signal corresponding to substantially 
the noise source alone; and 

generating an antinoise signal opposite in phase to said input signal, said output actuator means generating 
30 said acoustic output field from said antinoise signal. 

32. An acoustic noise control system for replacing the effects of a noise source with a desired signal, comprising: 

input transducer means for sensing the acoustic noise field generated by the noise source and for generating 
35 an input signal therefrom; 

output actuator means for generating an acoustic output field that is effective to produce the level and the 
characteristic of a replacing signal field; 

correction means for adjusting the input signal generated by said input transducer to compensate for the non 
linear characteristics of said input transducer and output actuator; 
40 echo cancellation means for removing from the input signal a portion of the output of said output actuator 

means fed back through said input transducer means, the output of said echo cancellation means representing 
a signal corresponding to substantially the noise source by itself; 
storage means for storing data representing said desired signal; 

replacing means for producing a replacing signal that is effective to replace said input signal with said desired 
45 signal; and 

wherein said input transducer means is located in relatively close proximity to said output actuator means. 

33. An acoustic noise control system for replacing the effects of a noise source with a desired signal, comprising: 

so input transducer means for sensing the acoustic noise field generated by the noise source and for generating 

an input signal therefrom; 

output actuator means for generating an acoustic output field that is effective to produce the level and the 
characteristic of a replacing signal field; 

correction means for adjusting the input signal generated by said input transducer to compensate for the non 
55 linear characteristics of said input transducer; 

input decoding means for removing extraneous signals from said input signal so as to generate a signal cor- 
responding to substantially the noise source alone; 
storage means for storing data representing said desired signal; 



21 



10 



15 



20 



25 



30 



35 



40 



45 



50 



55 



EP 1 414 021 A1 

replacing means for producing a replacing signal that is effective to replace said input signal with said desired 
signal; and 

wherein said input transducer means is located in relatively close proximity to said output actuator means. 
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FIG. 2 
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